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Preface 

Most  of  my  time  during  the  past  six  months  has  been  spent  preparing 
for,  working  at,  and  finally  completing  this  project.  As  time  progress- 
ed, I managed  to  develop  at  least  a measure  of  proficiency  in  working 
with  the  six  degree  of  freedom  equations  of  motion,  was  introduced  to 
the  utility  of  aerodynamic  coefficients,  and  gained  considerable  exper- 
ience working  with  a computer  simulation.  This  learning  experience,  un- 
like most  of  my  rote  textbook  encounters,  was  found  to  be  satisfying,  re- 
warding, and  real.  I also  found  the  professional  people  I came  in  con- 
tact with  during  the  project  to  be  knowledgable , motivated  engineers 
willing  to  share  their  time  and  talent.  My  association  with  my  advisor, 
Iir.  D.W.  Breuer  of  the  Aeronautics  and  Astronautics  department,  Air  Force 
Institute  of  Technology;  Major  Ed  Mirmak,  Air  Force  Avionics  Laboratory 
who  proposed  the  project;  and  Captain  Bill  Miklos,  Flight  Dynamics  Lab- 
oratory who  provided  the  computer  program,  made  the  project  not  only  pos- 
sible but  enjoyable.  An  enormous  amount  of  expertise  was  available  to  me 
through  these  people  as  well  as  from  those  who  have  previously  documented 
their  efforts  in  this  field.  This  project  represents  the  conclusion  to  a 
personal  challenge  of  the  first  degree  with  its  own  reward  at  the  time  of 
completion.  If  any  part  of  it  is  ever  useful  to  anyone  else,  then  all 
the  effort  is  even  more  worthwhile. 
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ABSTRACT 

A six  degree,  of  freedom  computer  simulation  was  used  to  Investigate 
the  lateral  progression  of  a free-fall  ballistic  trajectory  due  to  spin 
rate,  Magnus  aerodynamic  coefficients  and  initial  projectile  pitching 
motion.  The  increased  spin  rate  extends  the  projectile  impact  point 
both  downrange  and  cross  range  due  to  a slight  increase  in  time  of  flight 

i 

generated  by  a predominately  positive  angle  of  attack.  The  Magnus  force, 

Magnus  moment  and  side  force  coefficients, under  normal  release  conditions, 
presents  only  a minor  influence  and  can  be  omitted  from  the  simulation 
without  altering  the  trajectory  appreciably.  If  the  release  altitude  is 
sufficiently  high,  however,  the  small  influence  of  these  coefficients 

I 

could  propagate  to  a correctable  magnitude.  Projectile  oscillations  en- 
countered at  release  can  increase  the  lateral  progression  of  the  trajec- 
tory significantly,  especially  for  high  speed  deliveries.  Oscillations 
induced  by  an  initial  pitch  rate  of  the  projectile  generates  consider- 
ably more  lateral  deviation  than  an  initial  pitch  displacement  of  equiv- 
alent  maximum  amplitude.  Ibis  increased  lateral  displacement  causes  a 
corresponding  decrease  in  range  impact.  Other  observations  concerning 
the  performance  of  a projectile  following  a ballistic  trajectory  are  in- 
cluded as  support  material  in  the  last  section  of  the  study. 
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MAGNUS  EFFECTS  ON 
BALLISTIC  TRAJECTORIES 

I.  Introduction 

Purpose 

This  study  uses  computer  simulations  to  Investigate  three  factors 
that  effect  a free-fall  ballistic  trajectory.  The  first  objective  is 
to  determine  the  relative  Influence  of  Magnus  force , Magnus  moment,  and 
side  force  coefficients.  The  next  area  of  Interest  is  the  effect  of  os- 
cillatory motion  Induced  by  initial  pitch  rate  motion  compared  to  ini- 
tial pitch  displacement  of  the  projectile.  The  third  set  of  simulations 
determines  how  spin  rate  induced  by  the  projectile  fins  alters  the  tra- 
jectory. 

The  development  of  the  equations  of  motion  used  in  the  computer 
program  is  included  in  this  study,  followed  by  a discussion  of  the  aero- 
dynamic coefficients.  The  results  of  the  simulations  are  presented  in 
the  last  section. 

Background 

A ballistic  projectile  is  designed  to  spin  so  that  a predictable  tra 
jectory  is  generated.  Without  spin,  aerodynamic  lift  created  by  non-sym- 
metrical  shape  resulting  from  manufacturing  imperfections  would  cause  the 
projectile  to  fly  away  from  the  expected  ballistic  trajectory.  Spin  is 
also  required  for  stability  considerations.  If  the  projectile  has  in- 
sufficient angular  momentum,  it  will  eventually  begin  processing  and,  if 
the  spin  rate  is  in  the  neighborhood  of  the  natural  pitch  frequency,  it 


1 


—— 


HMPvnmQ! 

— 


AFIT/GA/ AA-77D-8 

■ay  develop  catastrophic  yaw.  Conversely,  with  excessive  angular  momen- 
tum, the  projectile  spin  axis  tends  to  become  gyroscopically  oriented 
with  respect  to  Inertial  space  and  will  not  track  properly  throughout 
the  trajectory.  Conventional  bombs  are  therefore  designed  so  that  a 
spin  rate  is  achieved  that  will  maintain  satisfactory  stability  within 
either  limit  of  spin  rate  (Ref  4i64  & 12i29).  This  spin,  while  desir- 
able and  necessary  for  a predictable  trajectory,  induces  lateral  errors 
which  are  attributable  to  Magnus  lift  and  a gyroscopically  induced  la- 
teral angle  of  attack.  These  effects  are  additive  and  can  cause  the  bomb 
to  progressively  deviate  laterally. 

Assumptions 

The  set  of  equations  that  completely  describe  the  trajectory  of  a 
spinning  ballistic  projectile  are  fully  coupled,  non-linear  expressions 
that  iu:e  developed  in  dynamics  texts  (such  as  Ref  4»99)»  The  solutions 
to  this  set  of  equations  requires  extensive  digital  computer  capability 
even  with  the  following  simplifying  assumptions! 

1.  The  earth  curvature,  variation  in  terrain,  and  earth  rotation 
are  neglected.  For  normal  conventional  weapon  delivery  airspeeds  and 
altitudes,  no  significant  error  is  introduced  since  the  range  and  time 
of  flight  of  the  projectile  is  sufficiently  small.  This  assumption  eli- 
minates the  earth  reference  and  allows  the  advantage  of  direct  transfor- 
mations between  the  body  and  inertial  reference  frames  (See  Fig  l). 

This  greatly  decreases  the  required  computer  time  with  no  appreciable 
loss  in  accuracy. 

2.  The  gravitational  field  is  considered  uniform  over  the  entire 
trajectory. 

3.  The  coriolls  acceleration  is  neglected  even  though  it  can  have 
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* small  but  disceznable  effect  on  the  projectile  Impact  point.  Fox  this 
investigation,  the  change  In  trajectory  Is  of  primary  Interest,  not  the 
precise  trajectory  Itself.  Coriolis  acceleration  Is  addressed  In  Ap- 
pendix E but  Is  not  written  Into  the  computer  program. 

4.  A standard  atmosphere  is  assumed  adequate  for  normal  release 
altitudes.  The  ARDC  1959  Model  Atmosphere  was  used  and  no  winds  were 
Included  in  the  simulation. 
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II.  Mathematical  Theory 


Reference  Frame  Transformation  Matrix 


The  orthogonal  reference  frame  used  In  this  development  Is  the  body 
fixed  reference  frame  (XYZ)  moving  with  respect  to  the  inertial  frame 
(xys).  The  origin  of  the  body  frame  is  located  at  the  center  of  gravi- 
ty of  the  projectile. 


Fig.  1.  Reference  Frames 


To  facilitate  the  direction  cosine  matrix  development,  the  body  ref 
erence  frame  is  rotated  90°  In  the  negative  direction  about  the  X axis 
to  allign  the  YZ  axes  with  the  corresponding  yz  axes.  The  following  se- 
quence of  Euler  angles  is  then  used  to  generate  the  direction  cosine  ma- 
trix (Ref  2i8)t  yaw  rotation  Y about  the  y axis,  pitch  rotation  6 about 
the  new  s'  axis,  then  roll  rotation  6 about  the  new  x"  axes. 
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The  resulting  matrix  ist 


CV 

0 

sr 

0 

1 

0 

-SY 

0 

CV' 

10  0 C9  SB  O 

0 C0  S<t)  -90  CQ  0 

0 - 50  C 0 0 O I 


which  expands  to  a single  Euler  angle  transformation  so  that  the  set  of 
equations  describing  the  velocity  vector  is  (Ref  7»10)i 


* cr  ca  so  co  sr  U 

y - -sesY-c^crse  c0ce  sficv  - cd>  svse  v 

£ -cdsr+s(t>cvse  -s&co  cficv  + sfisyso  w 


Where  x,  y,  z are  the  velocity  components  in  the  inertial  reference 
frame  and  u,  v,  w are  the  velocities  along  the  body  frame  (Pig  2).  S,  C 
denotes  sine  and  cosine  respectively. 

Equations  of  Motion 

The  relationship  between  vector  derivatives  in  any  two  reference 
frames  (where  Z denotes  any  vector)  can  be  expressed  as  (Ref  6;109)i 


-i-t  _ 

Z = Z + u x Z 


Then,  for  a velocity  vector  V 


jJL  ^_b  _b£  _ 

V = V + CJ  X V 
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or*  equivalently* 


- rv” 
m ru  - fw 
■pv  - <?  u 


where  the  components  of  F sure  expressed  as  the  combination  of  aerodyna 
mic  and  weight  forces. 


2 

where  Q = i = dynamic  pressure 
S = cross  sectional  area 
Cx,  Cy,  Cz  = aerodynamic  force  coefficients 

This  set  of  equations  (Ref  2 til)  describe  the  three  translational  de- 
grees of  freedom  of  a projectile  as  it  proceeds  along  its  path.  Three 
additional  equations  are  needed  to  account  for  the  spin,  pitch  smd  yaw 
motion  that  will  be  encountered. 

The  expression  relating  the  derivative  of  the  angular  momentum 
vector  between  the  body  and  Inertial  frame,  from  Eq  2,  is 


.2-1  -ib  _b.'  _ 

h - V\  + (J  x h 
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The  rate  of  change  of  angular  momentum  expressed  In  the  body  frame , 
for  constant  mass  projectiles , is 


The  origin  of  the  body  frame  of  reference  Is  located  at  the  center 
of  gravity  of  the  projectile  which  eliminates  the  products  of  inertia 
since  the  body  frame  is  located  on  principal  axes  and  the  moment  of  In- 
ertia about  the  y & z axes  are  equal.  The  scalar  components  of  Eq  9 
can  therefore  be  written  as 

& = mx  = + ?r[lz  - IY] 

lnY  = My  = Jyi  + T<f[  I*  -Tz]  (12) 

K*  = Mz  = x^r  + f?[ly  -T„] 

Since  I s I for  a projectile  with  rotational  mass  symmetry,  this  set 
y z 

is  equivalently  presented  as  (Ref  2:ll)t 

\f~\  (m* 

^nL  [mjt 

where  the  M vector  is  the  moment  produced  by  aerodynamic  forces  express- 
ed as  a function  of  the  aerodynamic  moment  coefficients  C^»  G^, 

(Ref  2tl2)t 
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’m* 

My 

- QSd 

i 

J X 

0 u 

1  

• 

Ch  J 

(14) 


where  d Is  the  diameter  of  the  projectile. 

Eqs  7 and  13  completely  describe  the  motion  of  a six  degree-of- 
freedom  ballistic  projectile  and  form  the  basis  of  the  six  degree-of- 
freedom  computer  program  used  In  this  Investigation. 


Differential  Equations  of  Motion  In  Rotational  Parameters 

Euler  angles  provide  an  effective  description  of  body  orientation 
in  space  and  the  heading,  pitch  and  roll  angles  are  easily  visualized. 
The  transformation  matrix  of  nine  direction  cosine  elements,  however, 
requires  considerable  computer  time  to  process  the  numerous  trigono- 
metric functions.  Also,  singularities  in  the  differential  equations  a- 
rlse  when  the  pitch  angle  of  ± 90°  is  reached,  as  well  as  inducing  un- 
acceptable truncation  error  when  integrating  in  the  neighborhood  of  the 
singularity  point.  The  use  of  Euler  rotational  parameters  (sometimes 
referred  to  as  quaternions)  provide  a computational  device  that  avoids 
the  singularity  limitation  of  Euler  angles.  This  system  uses  four  pa- 
rameters to  fix  the  position  of  a body  in  space i three  direction  co- 
sine angles  to  specify  the  orientation  of  the  spin  axis,  and  another 
rotational  parameter  to  specify  the  amount  of  spin  about  that  axis.  In- 
stead of  working  with  nine  direction  cosine  elements  and  six  constraint 

equations  from  the  Euler  set,  the  problem  is  now  reduced  to  expressions 

• • • 

for  the  rate  of  change  of  the  four  quaternion  parameters  ( A,,  Ax,  Aj, t 

4 ^223 

A0)  and  a single  constraint  equation  A,  + Aa  ♦ Aj  * A„  = I . The  com- 
puter time  required  to  solve  this  set  lo  reduced  but  still  remains  ex- 
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tensive.  The  rotational  parameter  equivalent  of  the  direction  cosine 
matrix  of  equation  2 Is  (Ref  2<7)i 

2(  I 2(A.Aj"*AoAa) 

2 ( A,  A,  ~ A,  A,)  2(A*+A,)-J  2(\\+\\) 

2 ( A,A3+  \\)  2 ( AaA5-  A.A, ) 2(  Ai  '+  A^)  - I 


(15) 


where 


\--ClC|C^ 

+ 

S-|5f 

A,  - C-f  S-f  S-£ 

S-£  C*  C.J 

At  = C-f  CIS^ 

+ 

Sf  s_&  cj£ 

(16) 

A3=C4  S ® G£ 

SlC|S{ 

The  rate  of  change  of  these  parameters  are  given  "by  (Ref  2tl0)i 


2 


O 

P 

? 


-r  - 


P 

o 

r 

? 


<?  r 

~r  ? 

o -p 

P o 


A o 

A, 

Aa 

(i?) 


The  algebra  required  to  develop  equations  15  and  16  is  quite  extensive 
as  shown  In  Appendix  D.  The  derivation  of  equation  17  Is  found  in  Ap- 
pendix B of  reference  2,  page  11. 

The  differential  equations  of  motion  are  now  in  the  programmed  form 
and  are  summarized  in  Table  I. 


t 

t 
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Table  I 

Differential  Equations  of  Motion 


The  kinematic  relationships  are i 


n 

X 

V 

— 

• 

'i 

•J 

■ 

** 

•< 
i 

- i 

2(  A,A2 - V^> 


2(  A,A,+  \A3) 

2(  A^  + A2)  - / 

2(Vh“  *<>*,) 


2(  A,  A,-  A.AJ 

2(  ^ + A^)-| 


X 


o 

-p 

-r 


P 

o 

r 


1 r 
-r  <? 


o 

P 


-P 

o 


A* 

A, 

A* 

A, 


The  dynamic  relationships  ares 


r . - 

u 

Fx 

^ w - r v 

m 

• 

V 

— 

fy 

- m 

ru  - pw 

• 

V 

a « 

A 

pv  - 

• — 

1 

r* 

PxP 

M 

*y7 


M- 


-dy-^)p 


- r 


where  F and  M are  described  In  Eqs  8 and  14. 
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III.  Method  of  Problem  Solution 


The  Computer  Program 

Solution  of  the  set  of  coupled,  non-linear  differential  equations 
listed  in  Table  I requires  extensive  computer  capability.  The  program 
used  in  this  study  was  originally  coded  for  the  Naval  Ordinance  Research 
Computer  and  has  been  updated  as  more  advanced  generations  of  computers 
became  operational.  The  Fortran  IV  program  used  for  this  study  was  de- 
veloped by  Charles  W.  Ingram  and  R.S.  Eikenberry  using  Cohen  & Werners' 
work  (Ref  2)  as  a point  of  departure.  It  was  written  originally  for 
stability  analysis  applications  rather  than  precision  trajectory  pre- 
diction. The  six  degree  of  freedom  computer  program  numerically  inte- 
grates the  set  of  kinematic  and  dynamic  Euler  equations  of  Table  I in 

I 

the  body  reference  frame  and  outputs  in  the  inertial  reference  frame. 

The  accuracy  of  the  program  is  limited  primarily  by  the  quality  of  the 
aerodynamic  coefficients  that  are  input.  The  complete  program  is  list- 
ed in  Appendix  F. 

Both  the  Runge-Kutta  numerical  integration  step  size  (value  of  A- 
(20)  Table  B-2)  and  the  output  time  increment  can  be  adjusted  by  the 
user  as  a means  of  tailoring  the  computer  time  and  cost  against  the  ac- 
curacy required.  If  high  projectile  spin  rates  are  encountered,  coarse 
step  sizes  on  the  order  of  .01  second  could  result  in  erroneous  outputs 
due  to  insufficient  data  points  to  properly  resolve  the  integration 
problem.  The  time  step  should  be  small  enough  to  provide  at  least  ten 
integration  steps  during  one  revolution  of  the  highest  spin  rate  expect- 
ed. A step  size  of  .002  second  is  used  for  most  simulations  in  this  re- 
port. Since  extensive  core  memory  and  computer  time  is  required  to  pro- 


a 
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cess  this  program , it  is  written  for  only  a single  trajectory  per  run. 

A core  nemory  of  100,000  was  adequate  for  all  trajectories  used  in  this 
study.  (An  integration  step  size  of  .001  second  requires  in  excess  of 
200  seconds  of  computer  time  for  a 10,000  ft  drop.) 

Discussion  of  the  computer  program  input  and  output  is  included  in 
Appendices  B and  C. 

Projectile  Model  Description 

The  standard  Mark  82  warhead  with  an  experimental  tail  section  was 
the  projectile  model  used  to  determine  the  aerodynamic  coefficients. 

Ihe  cylindrical  aft  section  has  four  slotted  fins,  each  fitted  with  a 
trailing  edge  wedge  to  provide  the  roll  driving  moment.  The  designation 
of  the  configuration  is  FFSW  (fin,  fixed,  slotted,  wedge).  This  particu- 
lar projectile  was  selected  because  of  the  similarity  to  existing  conven- 
tional weapons  and,  significantly,  most  of  the  aerodynamic  coefficients 
were  available  (Ref  14). 
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PIP 


- 


The  cylindrical  tail  section  provides  both  a higher  drag  coeffi- 
cient and  a housing  for  an  inflatable  high  drag  device  that  may  be  in- 
stalled; the  slotted  fins  increase  the  probability  of  a well-behaved 
roll  rate  throughout  the  trajectory  (Ref  3»1). 

A wedge  angle  of  ten  degrees  is  used  for  most  of  the  simulations 
in  this  study.  This  wedge  angle  provides  sufficient  roll  acceleration 
to  overcome  adverse  dynamics  encountered  during  release.  The  steady 
state  spin  rate  produced  is  fast  enough  to  Investigate  the  influence 
of  Magnus  coefficients  without  inducing  the  type  of  instabilities  that 
can  occur  when  the  spin  rate  exceeds  approximately  ten  times  the  nuta- 
tion frequency  (Ref  8;  Sec  3,5). 

Aerodynamic  Coefficients 

Values  of  the  aerodynamic  coefficients  used  in  equations  8 and  14 
must  be  adequately  determined  by  wind  tunnel  testing  of  the  projectile 
model  before  satisfactory  trajectory  simulations  are  possible.  The 
aerodynamic  force  and  moment  coefficients  can,  for  most  applications, 
be  reduced  to  the  expressions  listed  in  Table  II.  These  expressions  in- 
clude cross- velocity  terms  in  both  the  force  and  moment  coefficients, 
but  cross  spin  terms  are  included  only  in  moment  coefficients  about 
pitch  and  yaw  axes.  (See  Fig  6 for  definitions  of  cross- velocity  and 
cross  spin  terms.)  Table  III  describes  the  terms  which  make  up  the  total 
aerodynamic  force  and  moment  coefficients.  These  mathematical  expres- 
sions were  determined  in  Cohen  and  Werners'  work  (Ref  2il5-19). 

In  order  to  illustrate  the  utility  of  the  aerodynamic  coefficient 
terms,  the  expression  for  Induced  roll  moment  is  used  as  an  example 
(See  Table  III)j 
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Cj0  ♦ C I?  sin  rjf*  + Cjg  cos  rj 


where 


* Angle  to  cross  velocity  vector  (See  Fig  6) 

C|0  = Initial  curve  shift  (See  Fig  5) 

Cj| 7 • Cjg  = Amplitude  of  curve 
rj  = Number  of  fins 

The  CjJq  term  is  the  residual  rolling  moment  that  remains  after  the  pro- 
jectile has  been  aligned  in  the  wind  tunnel  at  the  angle  of  attack  of 
interest  with  no  fin  deflection j that  is,  an  index  of  rolling  moment 
due  to  body  shape. 

The  next  set  of  terms,  Cj?  sin  r)  0'  + Cjjg  cos  q ' , is  the  mathema- 
tical format  required  to  describe  the  data  curve  obtained  from  the  wind 
tunnel  measurements.  This  application  can  best  be  appreciated  by  the 
illustration  of  a spinning  projectile  at  various  stages  of  roll  moving 
away  at  an  angle  of  attack  sufficient  to  cause  a vortex  to  trail  along 
the  top  of  the  missile  as  in  Figure  4. 


t • 22.5 

• • *5 

• • W.5 

Positive  Roll 

■eg.tlYt  Roll 

ZnflooAco 

Xaflucoc. 

Fig.  4.  Periodic  Change  in  Aerodynamic  Coefficient  Cj 


At  the  beginning  of  the  sequence,  the  trailing  vortex  is  cut  symmetri- 
cally by  the  fin,  so  no  net  roll  influence  is  imposed  by  the  dynamic 
pressure  change  on  the  fin  due  to  the  vortex.  After  22.5  degrees  of 
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roll,  the  vortex  Is  in  a position  to  exert  a maximum  rolling  influence. 
At  the  45°  angle,  a condition  of  symmetry  is  again  encountered.  For 
this  type  of  oscillating  characteristic,  Cj  could  be  graphed  for  the 
entire  sequence  as  the  solid  curve  of  Fig  5.  The  graph  obtained  from 
the  actual  wind  tunnel  data  for  the  coefficient  of  Interest,  however, 
is  generally  skewed  from  the  theoretical  curve  as  Indicated  by  the  dot- 
ted curve. 


Fig.  5.  Theoretical  vs  Experimental  Data  Curves 


The  three  terms,  C^Q  ♦ C slnr}*'  ♦ C|Q  costr^',  are  therefore 
used  in  whatever  combination  of  magnitudes  and  sinusoidals  that  produce 
the  best  simulation  of  the  actual  wind  tunnel  data.  The  relative  mag- 
nitude of  Cjpo  versus  the  amplitude  of  the  sinusoidals  often  makes  one 
or  the  other  of  the  terms  negligible. 

The  remaining  trigonometric  expressions  in  Table  II  not  defined  in 
Table  III  are  attributable  to  transformations  between  reference  frames. 

Most  of  the  aerodynamic  coefficient  arrays  used  in  this  computer 
program  had  been  previously  determined  through  a joint  effort  between 
the  Naval  Surface  Weapons  Center,  Dahlgren  Laboratory,  and  the  Aero- 
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where  (Ref  2) 


q*  - 

^2  + r2 

, magnitude  of  cross-spin  (rad/sec) 

a *» 

arc  tan , 

I 4*  v/2 

/ , magnitude  of  yaw  or  angle  of  attack  (deg) 

v . u. 

1 

4?  « 

arc  tan 

V 

— - (angle  about  X-axis  from  cross-velocity  vector 

1 

of  llie  center  of  gravity  of  missile  to  Z fir}  (deg) 


arc  tan  angle  about  the  X-axis  from  cross-spin  vector 

to  Y fin)  (deg) 


1? 


I HM 


Force  And  Moment  Coefficients  (Adapted  from  Ref  2) 
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TABLE  III 


Aerodynamic  Coefficients 
(Adapted  from  Ref  2) 


Description  of 

Coefficient 

Array 
Symbol  # 

CX 

St 

Axial  force 

Cy0  ♦ Cy7s!nij0*  +Cygcosij£' 

Induced  side  force 

CY7 

S 

Magnus  force 

CYP 

CjjQ+C^y  sinn^'  +CZ8cos^’ 

Normal  force 

CZ 

CyW 

Trim  force  along  Y-axls 

N/A 

<*(«) 

Trim  force  along  Z-axis 

N/A 

Cj^+C^  s In n<t>'  + Cigcosij^* 

Induceu  roll  moment 

CLGA 

Cj(iw) 

roll  moment  due  to  fin  cant 

CLDW 

% 

Roll  damping  moment 

CLP 

Cm0+Cm■]sinrl4,,  ^CmgCos^' 

Restoring  moment 

CM 

^"qo+Cmq7  sin^”+cmqgCOST,^> 

'0f£ 

a 

Damping  moment 

CM^ 

Cno+cn7  sln  ^,+cnj  cos'>*’ 

Induced  side  moment 

CNGA 

% 

Magnus  moment 

GRP 

Cm« 

Trim  moment  about  y-axis 

CMDS 

Ca(0 

Trim  moment  about  z-axis 

N/A 

4 <CJtftf0+Cif0.s,n<^' + C^«2  C0S^ 

Roll  moment  due  to  eccentric  tall 

N/A 

* See  Appendix  G 
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space  Research  Laboratory , Wright-Patterson  Air  Force  Base  and  were 
available  for  this  project  (See  Appendix  G).  The  Magnus  force  coeffi- 
cient, however,  was  not  available  and  had  to  be  determined  so  that  its 
effect  on  lateral  error  could  be  studied.  The  Magnus  force  coefficient 
is  determined  from  the  definition 


8c 


8 (pd/2V) 


(23) 


I 


where  p = spin  rate,  d = projectile  diameter,  V = wind  tunnel  velocity. 

This  relation  states  that  the  slope  of  a line  fit  through  a plot 
of  Cy  versus  pd/2V  should  yield  the  nominal  value  of  Cy^  for  a specific 
angle  of  attack  and  Mach  number.  Tables  A-I  and  A-II  show  the  wind  tunnel 
data  taken  at  Mach  0.8  at  an  angle  of  attack  of  1.8?  degrees.  Cy  was 
plotted  against  pd/2V  and  then  a "best  fit"  line  was  drawn  through  the 
data  points  (Fig  8).  Greater  weight  was  given  to  the  lower  values  pd/2V 
so  that  the  cluster  of  data  points  in  the  neighborhood  of  the  steady 
state  spin  condition  would  not  overly  influence  the  slope.  No  attempt 
was  made  to  force  the  line  through  the  origin.  As  long  as  the  axis  in- 
tercept remains  small,  the  error  in  the  wind  tunnel  measuring  instru- 
ments will  dominate.  The  primary  objective  is  to  provide  a best  esti- 
mate of  the  most  representative  slope  of  the  plot.  The  slope  of  the 
resulting  line  produced  the  nominal  value  of  S that  was  entered  as  a 
single  element  in  the  CYP  array  corresponding  to  Mach  = 0.8  and  c*.  = 2 
(Ref  Table  IV).  The  remaining  ?1  elements  of  the  8x9  matrix  were  de- 
termined in  the  same  manner.  Note  that  the  Mach  = 0.4  values  are  used 
to  fill  the  Mach  r 0 row.  This  is  done  because,  at  lower  Mach  numbers, 
the  compressibility  effects  diminish  and  the  coefficient  remains  nearly 
constant.  Fig.?  shows  a typical  force  coefficient  progression. 
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0 .5  1.0  l.s 

MACU  NUMBER. 


7.  Force  Coefficient  vs  Mach  Number 


The  alpha  = 50  column  contains  unreliable  data  but  is  used  only  to 
provide  a value  for  the  program  to  interpolate  against  in  the  event  ©<« 
25°  is  exceeded.  The  computer  program  is  written  so  that  interpolation 
is  accomplished  for  any  intermediate  value  of  angle  of  attack  and  Mach. 
The  angles  of  attack  used  in  this  study  were  well  below  25  degrees  so 
this  ill-defined  edge  of  the  array  was  not  encountered. 


Approach  to  the  Investigation 

The  primary  objective  of  this  study  is  to  Investigate  the  relative 

contributions  of  the  Magnus  coefficients  to  lateral  deflection  of  a spin 

ning  projectile  following  a ballistic  trajectory.  The  general  approach 

is  to  isolate  the  effect  of  the  Magnus  force  coefficient  Cv  , the  in- 

yP 

duced  side  coefficient  Cv  , and  the  Magnus  moment  coefficient  Cn  on 

J7  P 

the  trajectories  initiated  from  identical  release  conditions.  Each  co- 
efficient is  studied  individually,  and  in  pairs;  then  compared  to  the 
performance  of  all  three  coefficients  together.  The  initial  conditions 
Include  a range  of  airspeed,  altitudes  and  attitude  pertubations . 

The  Magnus  force  coefficient,  Cv_ , is  of  interest  because  it  is 
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the  source  of  Magnus  lift  that  is  generated  by  a spinning  body  subject 


to  a relative  wind  component  acting  perpendicular  to  the  body  spin  axis 
This  force  will  therefore  be  generated  anytime  the  spinning  projectile 


is  falling  with  an  angle  of  attack.  The  direction  of  the  lateral  Mag- 


nus lifting  force  is  determined  by  the  direction  of  spin  and  the  cross- 


velocity component  w.  The  spin  rate,  p,  causes  a velocity  differential 


to  be  produced  as  shown  in  Pig.  9*  A corresponding  pressure  differen- 
tial is  generated  which  produces  lift  in  the  lateral  direction. 


contributes  to  the  deflection  of 


a projectile  by  influencing  side  lift  due  to  a lateral  angle  of  attack 


Pig.  9*  Magnus  Lift 


The  Magnus  moment  coefficient,  On  , can  directly  effect  lateral 

P 

progression  of  a projectile  following  a ballistic  trajectory.  Magnus 
moment  is  generated  when  the  projectile,  at  a sufficiently  high  angle 
of  attack,  trails  a wake  vortex  that  blanks  out  the  rolling  lift  of  the 


fin  passing  through  it.  This  causes  a change  in  roll  torque  and  a fin 
force  imbalance  that  produces  a yaw  moment  on  the  projectile  (Pig  10). 


The  value  of  the  aerodynamic  coefficient  is  a measure  of  the  sensiti- 


vity of  the  projectile  to  this  effect 
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Fig.  10.  Ma&ius  Moment 


The  release  velocities  used  as  initial  conditions  are  600  ft/sec 
for  a lower  limit  (approximately  350  knots)  and  1000  ft/sec  for  the  high 
speed  delivery.  The  lower  limit  was  selected  as  a representative  velo- 
city for  a slow  delivery.  The  1000  ft/sec  velocity  remains  in  the  sub- 
sonic region  and  provides  adequate  velocity  spread  for  comparisons.  A 
velocity  of  1200  ft/sec  was  used,  in  some  cases,  to  investigate  the  ef- 
fect of  a supersonic  delivery  (Mach  1.1  at  10,000  ft). 

The  10,000  foot  altitude  was  selected  for  most  of  the  simulations. 
This  altitude  provides  sufficient  time  of  flight  to  identify  any  signi- 
ficant trajectory  trends  and  does  not  require  prohibitive  computer  time 
to  process.  An  altitude  of  20,000  feet  was  used  when  a lower  air  den- 
sity trajectory  was  desired  for  comparison. 

All  releases  are  from  a level  flight  path.  The  projectile  is, 
however,  perturbed  initially  under  three  conditions!  zero  perturbation, 
a nose-down  rate  motion,  and  a nose-down  displacement.  The  zero  per- 
turbation condition  simulates  a projectile  perfectly  allied  with  the 
velocity  vector  at  release.  The  pitch  rate  motion  of  q » -.5  radians/ 
sec  is  used  to  investigate  the  effect  of  motion  that  could  be  induced 
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by  the  bomb  ejection  mechanism.  The  Initial  angle  of  attack  of  o^s 
-7.75  degrees,  the  third  condition,  is  approximately  the  same  maximum 
amplitude  achieved  by  the  preceding  initial  motion,  but  obtained  in  a 
different  manner. 

Spin  rate  of  the  projectile  directly  influences  Magnus  lift  and, 
consequently,  lateral  progression  of  the  trajectory.  Wedge  angles  of 
5,  10,  and  15  degrees  were  used  in  the  simulations  to  investigate  the 
relative  influence  of  the  different  wedges. 

The  parameters  studied  under  various  initial  conditions  axe  listed 
an  an  overview  in  Table  V. 
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IV.  Results  of  the  Simulations 

Effect  of  Magnus  Coefficients 

The  first  simulations  were  run  to  determine  the  relative  influence 
of  the  Magnus  coefficients.  Initial  release  conditions  of  level  flight 
at  an  altitude  of  10,000  feet  and  600  feet  per  second,  using  all  twelve 
available  aerodynamic  coefficient  arrays  produced  the  result  shown  in 
Fig.  11.  The  ground  zero  impact  point  for  the  trajectory  using  all  co- 
efficients was  5?. 73  feet  to  the  left  of  the  inertial  x-y  plane  contain- 
ing the  release  point.  With  three  coefficients  removed  that  influence 
lateral  displacement  (C^,  Cy^,  Cy^),  the  impact  point  was  changed  to 
56.58  feet.  This  total  spread  of  only  1.15  feet,  generated  over  a re- 
latively long  trajectory,  is  lost  in  the  accuracy  limitations  of  the 
aerodynamic  coefficients  and  the  simplifying  assumptions  used  in  the 
equation  of  motion  development.  This  small,  Magnus  induced  lateral  mo- 
tion is  initiated  early  and  propagates  slowly  throughout  the  trajectory. 

A second  set  of  simulations  were  run  at  a velocity  of  1200  ft/sec 
to  determine  if  higher  velocity  would  cause  these  coefficients  to  gen- 
erate a more  significant  change  in  the  trajectories.  The  net  difference 
in  the  trajectories  with  and  without  the  three  coefficients  remains  on 
the  order  of  only  a few  feet. 

There  are  a number  of  possible  explanations  why  the  Magnus  influ- 
ence is  so  small.  Magnus  lift  depends  on  both  the  relative  velocity 
component  acting  on  the  projectile  due  to  angle  of  attack  and  the  spin 
rate  of  the  projectile.  For  normal  weapon  deliveries,  the  angle  of  at- 
tack is  relatively  small  on  release  and  dampens  toward  zero  early  in  the 
fall.  During  the  time  the  angle  of  attack  is  at  a maximum,  shortly  af- 
ter release,  high  spin  rate  has  not  yet  developed.  After  the  projectile 
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am  ij.i 


Fig.  11.  Effect  of  Magnus  k Side  Force  Coefficients  on  Lateral  Displacement 
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has  fallen  long  enough  to  generate  a high  spin  rate,  the  angle  of  attack 
Is  negligible.  The  Magnus  side  force  due  to  body  lift  is  negligible  be- 
cause the  lateral  angle  of  attack  is  small.  The  Magnus  moment  is  also 
small  since  there  is  no  significant  wake  generated  at  small  angles  of 
attack  to  produce  a lift  differential  in  the  fins. 

If  oscillatory  motion  is  encountered  during  release,  the  projectile 
angle  of  attack  alternates  from  positive  to  negative.  The  Magnus  lift 
will  therefore  alternate  from  positive  to  negative  and  average  the  net 
influence  to  zero. 

Effect  of  Initial  Pitch  Motion 

To  simulate  an  initial  release  condition  more  realistic  than  per- 
fect alignment  of  the  projectile  with  the  velocity  vector,  an  initial 
angle  of  attack  of  -?.75  degrees  was  used  as  an  perturbation.  The  same 
release  initial  conditions  were  used  as  before  and  simulations  were  run 
first  using  all  coefficients,  and  then  with  Cy^,  C^,  Gy^  removed.  There 
was  essentially  no  difference  in  the  resulting  trajectories  for  o<q=  -7.75 
degrees  compared  to  c<o=  0 as  shown  in  Table  VI). 

An  initial  pitch  rotation  of  the  projectile  at  release,  as  might  be 
induced  by  an  ejector  mechanism,  was  then  implemented  by  entering  a val- 
ue of  qj,  = -.5  rad/sec.  This  initial  rotation  produced  the  same  maximum 
amplitude  of  oscillation  as  the  c<0=  -7*75°  case,  but  a considerable  dif- 
ference in  the  resulting  impact  point  of  the  projectile  was  generated. 

The  increased  lateral  displacement  caused  a corresponding  decrease  in 
range  which  is  consistent  with  the  principle  of  conservation  of  energy. 

The  oscillation  induced  by  an  initial  angular  displacement  produced  a 
trajectory  that  impacted  at  14,165  feet  downrange  and  57*73  feet  cross 
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Fig.  12.  Effect  of  Pitch  Motion  and  Magnus 
Coefficients  on  Lateral  Displacement 
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range  to  the  left.  The  rate  Induced  oscillation,  with  the  sane  release 
conditions,  inpacted  at  14,1?9  feet  downrange  and  73.02  feet  cross  range 
(See  Table  VI  and  Fig  12). 

tb  Investigate  the  cause  of  the  difference  in  trajectory,  the  nose 
track,  viewing  the  projectile  from  behind,  was  plotted  for  the  two  dif- 
ferent Inputs  in  Fig.  13.  The  beta  versus  alpla  plot  represents  a close 
approximation  of  the  nose  track  of  the  projectile  early  in  the  trajec- 
tory before  the  body  reference  frame  and  the  projectile  velocity  vector 
rotate  a significant  amount  with  respect  to  the  inertial  frame.  The  mo- 

a 

tion  of  the  projectile  nose  for  the  two  conditions  is  quite  different. 

Fig  13a  represents  single  arm  coning  motion  caused  by  nutation.  No  pre- 
cession is  evident  and  the  damping  is  relatively  constant.  The  nose 
track  in  Fig  13b  shows  less  coning  and  displays  two  arm  motion  with  both 
nutation  and  precession.  This  oscillation  dampens  faster  in  pitch  than 
in  yaw.  The  direction  of  nutation  is  opposite  to  the  projectile  spin 
rate  while  the  precession  is  in  the  same  direction.  Figs  13c,  d,  e,  f 
were  plotted  to  identify  any  parameters  that  could  cause  the  difference 
in  trajectories.  These  plots  show  that  the  difference  in  spin  rate  build- 
up and  velocity  progression  is  very  small.  A slight  phase  shift  in  yaw 
oscillation  is  apparent,  and  the  initial  slopes  of  the  ^ versus  t plots 
are  noticably  different.  The  last  plot  shows  a significantly  different 
behavior  in  the  angle  of  attack  oscillations  between  the  two  initial  con- 
ditions . 

The  lateral  displacement  for  both  initial  conditions  is  plotted  a- 
galnst  time  in  Fig.  14  to  further  define  the  projectile  motion  shortly 
after  release.  The  linearized  average  slope  of  these  plots  produces  a 
measure  of  the  initial  lateral  average  velocity.  The  slope  for  the 
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Fig.  13.  Projectile  Nose  Track  After  Release 
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e<*-7.?5°  Initial  condition  is  approximately  1 ft/secj  for  q • -,j  rad/sec, 
the  slope  is  close  to  1.5  ft/sec.  These  values  of  slope  are  significantly 
different  and  can  be  modeled  as  an  initial  lateral  velocity  of  the  pro- 
jectile at  release.  The  Initial  lateral  velocity,  vq,  can  be  input  into 
the  simulation  to  determine  its  effect.  Fig.  15  shows  the  lateral  pro- 
gression of  the  trajectories  resulting  from  an  initial  side  velocity  of 
-5,  0,  +5  ft/sec.  The  remaining  initial  release  conditions  were  the 
same  for  all  three  trajectories.  As  the  plots  show,  the  magnitude  of  vq 
directly  influences  the  overall  trajectory.  An  equation  for  this  lateral 
initial  velocity  is  given  by  (Ref  li 53-68 )i 


where  C„  B the  normal  force  coefficient  (anting  in  plane  containing  the 
relative  wind). 

A = frontal  area  of  the  projectile, 
q s dynamic  pressure. 

H = mass 

T ■ time  constant  for  oscillation  damping. 

0J*  frequency  of  oscillation. 

Q s yaw  angle  of  attack. 

£ a yaw  rate. 

If  the  Influence  of  the  displacement  term  is  considered  small. 


with  respect  to  the  rate  term,  the  expression  reduces  to i 
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Fig.  15.  Effect  of  Initial  Side  Velocity  on  Lateral  Progression 
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1 

CM-  H d 


where  » moment  coefficient  (with  respect  to  total  angle  of  attack). 


I * transverse  moment  of  inertia. 


d x max  diameter  of  projectile. 


and  the  following  relation  has  been  usedt 


1* 

1 


°M  * <1  d 


Equation  25  shows  that  the  lateral  velocity  induced  by  oscillation 


of  the  projectile  is  sensitive  to  the  normal  force  coefficient  Cz,  the 
moment  coefficient  C^,  and  the  yaw  rate  <£> . 

The  magnitude  of  the  side  velocity  predicted  by  this  closed  form 


approximation  is  given  by  (Ref  1*65) 


V b yr-  (.362)  ft/sec/degree 

Z SjJ 


where  the  constant,  0.362,  is  the  result  of  equation  24  (displacement 
term  neglected)  using  a velocity  of  750  ft/sec,  an  altitude  of  $000  ft, 


and  coefficients  of  a projectile  similiar  to  that  in  this  investigation 


(coefficients  used  in  this  derivation  are  determined  with  respect  to 
the  total  angle  of  attack  instead  of  pd/2V).  These  same  initial  con- 


ditions were  used  in  a six  degree  of  freedom  simulation  and  the  lateral 


displacement  and  yaw  oscillations  were  plotted  in  Fig.  16.  From  these 


plots,  a representative  yaw  rate  and  frequency  of  oscillation  were 
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Fig.  16.  Lateral  Displacement  and  Yaw  Oscillations 
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estimated  as  O.349  radians/sec  and  5*3^  radians/sec  respectively.  The 
closed  form  approximation  (Eqn  2?)  then  produces  a value  of  ■ 1 .36 

feet  per  second. 

The  plot  of  lateral  displacement  versus  time  indicates  the  initial 
lateral  velocity,  according  to  the  simulation,  is  approximately  1.8  feet 
per  second.  Thi6  is  on  the  order  of  2 5%  difference  between  the  six  degree 
of  freedom  simulation  results  and  the  closed  form  expression.  This  dif- 
ference can  be  attributed  to  the  more  extensive  simplifying  assumptions 
used  to  make  the  closed  form  expression  possible. 

Effect  of  Spin  Rate 

The  effect  of  spin  rate  induced  by  various  wedge  angles  is  investi- 
gated at  different  airspeeds  and  release  perturbations  (Table  VIl).  The 
spin  rate,  induced  by  the  wedge  angle  of  fifteen  degrees,  produces  the 
most  nominally  linear  lateral  deviation  throughout  the  entire  trajectory. 
This  pattern  holds  for  both  velocities  investigated  (Fig  1?). 

Table  VII  also  reveals  that  increasing  the  spin  rate  of  the  projec- 
tile causes  the  trajectory  to  impact  longer  and  wider.  This  apparent 
contradiction  to  the  principle  of  conservation  of  energy  was  initially 
thought  to  be  an  error  caused  by  using  an  integration  step  size  that  was 
too  large.  A step  size  of  .002  second,  however,  is  adequate  for  the 
spin  rates  encountered. 

Another  possible  explanation  is  that  the  projectile  developed  gyro- 
scopic rigidity  and  was  starting  to  fly  at  a net  angle  of  attack  rather 
than  averaging  to  zero  angle  of  attack  in  its  oscillations.  This  con- 
cept is  supported  by  the  simulation  readout  which  shows  the  angle  of  at- 
tack of  the  projectile  with  a five  degree  wedge  angle  oscillates  on  the 


Effect  of  Spin  Rate 
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order  of  t 0.1  degree  in  the  latter  half  of  the  trajectory.  The  fifteen 
degree  wedge  readout,  however,  shows  only  about  half  the  magnitude  of 
oscillation  but  the  oscillations  are  predominantly  positive  in  the  last 
half  of  the  trajectory.  In  fact,  the  angle  of  attack  does  not  go  neg- 
ative at  all  during  the  final  ten  seconds  of  flight.  Also,  Table  VII 
Indicates  the  total  time  of  flight  of  the  projectile  spinning  at  the 
faster  rates  is  Increased  very  slightly  (on  the  order  of  0.02  second). 

The  downrange  velocity  just  prior  to  impact  is  approximately  750  ft/sec 
which  produces  an  opportunity  for  the  projectile  to  impact  15  feet  long- 
er in  the  extra  0.02  second.  This  is  compatible  with  the  data  shown  in 
Table  VII. 

Other  Observations 

Fig.  18  shows  the  aerodynamic  coefficients,  Cyy,  Cy^  A , have  no 
discemable  effect  on  the  projectile  velocity.  The  velocity  progression 
graphs  remained  the  same  regardless  of  whether  or  not  these  coefficients 
were  included  in  the  simulation. 

The  wedge  angle  did  not  change  the  velocity  progression  significantly. 
Simulations  modeling  spin  rates  as  5*  10,  or  15°  wedge  angle  indicates  the 
velocity  progression  is  not  influenced  by  spin  rate.  This  observation  is 
subject  to  the  projectile  spin  rate  remaining  low  enough  to  preclude  in- 
stability. 

The  higher  dynamic  pressure  encountered  during  release  at  high  air- 
speeds causes  a significantly  greater  velocity  change  of  the  projectile 
after  release.  The  magnitude  of  the  difference  in  the  release  and  steady 
state  velocity  of  the  projectile  governs  how  the  velocity  will  progress 
throughout  the  trajectory.  The  projectile  shows  significant  velocity 
change  throughout  the  trajectory  when  released  from  a low  speed  condition. 
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Fig.  18.  Effect  of  Release  Velocity  and  Altitude  on  Velocity  Progression 


AFIT/CA/AA-77D-8 


The  projectile  shows  ouch  less  acceleration  when  released  from  a high 
dynamic  pressure  condition. 

Figure  19  investigates  projectile  acceleration  from  release  to  im- 
pact. The  slope  of  the  velocity  vs  time  curve  shows  that  maximum  decel- 
eration is  proportional  to  the  drag  force  encountered  at  release  and,  as 
a result,  velocity  loss  is  greater  at  high  release  velocities.  The  time 
required  for  transition  ftom  deceleration  to  acceleration  also  increases 
with  higher  release  velocities.  The  same  plots  apply  for  each  of  the 
projectile  fin  wedge  angles  investigated.  Also,  there  are  no  signifi- 
cant changes  in  the  graphs  for  an  initial  angle  of  attack  of  o<0=  0 ver- 
sus c*=  -10. 

Fig.  20  shows  the  lateral  displacement  progression  of  a projectile 
subjected  to  two  different  release  perturbations  ( <X0=  -10°  and  qQ  = -.5 
rad/sec)  at  two  different  release  velocities.  The  slope  of  a pair  of  the 
Ottrves  is  nearly  identical  at  any  given  time  but  the  magnitude  of  the  la- 
teral displacement  attained  in  the  same  given  time  can  be  significantly 
different  at  the  lower  airspeed. 

Fig.  21  shows  that  the  rate  of  altitude  loss  does  not  change  signi- 
ficantly with  initial  release  airspeeds.  For  the  10,000  ft  drop,  time  of 
flight  through  a standard  atmosphere  is  increased  on  the  order  of  one 
second  compared  to  the  time  of  flight  through  a vacuum. 

The  effect  of  release  velocity  and  altitude  is  shown  as  a composite 
in  Fig.  22.  The  set  of  trajectories  initiated  at  10,000  ft  is  essential- 
ly a duplication  of  Fig.  11  and  is  included  with  the  20,000  ft  set  to 
give  an  Indication  of  the  effect  of  air  density  on  the  lateral  deflection. 
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Fig.  19*  Effect  of  Release  Velocity  on  Total 
Projectile  Acceleration 
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Fig.  21.  Effect  of  Release  Velocity  on  Altitude 
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V.  Conclusions 

1 

The  Magnus  force,  Marius  moment,  and  side  force  coefficients,  for  a 
projectile  of  the  type  modeled  in  the  simulations,  and  for  normal  release 
conditions,  exert  only  a very  small  influence  on  the  impact  point.  The 
cross  range  displacement  attributable  to  these  coefficients,  for  a tra- 
jectory initiated  at  10,000  feet,  are  on  the  order  of  only  a few  feet. 

The  influence  on  the  downrange  component  of  the  trajectory  is  also  very 
small.  These  coefficients  could  be  influential  under  conditions  of  high 
angle  of  attack  and  high  spin  rates  which,  for  normal  releases,  are  not 
present  simultaneously.  When  the  projectile  oscillation  and  the  corres- 
ponding angle  of  attack  is  maximum  following  release,  the  spin  rate  has 
not  yet  developed.  As  the  projectile  progresses  along  the  trajectory, 
the  spin  rate  builds  up  tut  the  oscillations  dampen  to  a very  small  value. 

Also,  the  oscillations  alternate  about  the  velocity  vector  of  the  center 
of  gravity  of  the  projectile  so  the  Magnus  force  direction,  being  depen- 
dent on  the  direction  of  the  angle  of  attack,  also  alternates  and  tends 
to  average  to  zero  effect. 

The  lateral  displacement  of  the  trajectory  is,  in  general,  increas- 
ed with  increasing  spin  rate.  Insufficient  spin  rate  can  be  a stability 
factor  as  can  excessive  spin  rate.  Low  values  of  angular  momentum  can  al- 
low the  projectile  to  develop  catastrophic  yaw.  High  spin  rates  induce 
gyroscopic  resistance  to  change  in  projectile  orientation,  allowing  a 
positive  angle  of  attack  to  build  up  as  the  trajectory  progresses.  A net 
lift  is  then  developed  which  increases  the  time  of  flight  of  the  projec- 
tile end,  consequently,  the  cross  range  and  downrange  impact  point. 

A projectile  released  at  zero  angle  of  attack  with  no  initial  per- 
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tur bat Ions  produces  the  least  lateral  displacement  of  the  trajectory.  If 
the  projectile  Is  released  with  an  Initial  negative  angle  of  attack,  the 
lateral  progression  Is  Increased  only  slightly.  If  the  projectile  Is  re- 
leased with  a rate  rotational  motion,  as  could  be  Induced  by  the  ejector 
mechanism  or  an  initial  positive  angle  of  attack,  the  lateral  progression 
is  significantly  Increased.  Dlls  Increase  is  directly  proportional  to 
the  amplitude  of  the  projectile  oscillations.  The  difference  In  lateral 
progression  between  the  rate  Induced  and  displacement  Induced  oscillation 
is  less  at  higher  airspeeds.  Also,  as  oscillations  Increase  the  cross 
range  impact  point  of  the  trajectory,  there  is,  in  general,  a correspond- 
ing decrease  in  the  downrange  impact  point.  These  results  suggest  the 
following  conclusions! 

1.  Magnus  force,  as  such,  exerts  no  significant  influence  on  the 
trajectory  of  a projectile  such  as  the  Mk-82  bomb.  For  most 
release  airspeeds  and  altitudes,  compensation  for  this  specific 
effect  in  a weapons  delivery  system  would  not  be  cost  effective. 

2.  The  projectile  should  be  designed  to  spin  as  slow  as  possible 
consistent  with  stability  considerations. 

3.  Release  mechanisms  should  be  designed  so  no  significant  moment 
is  applied.  If  assisted  separation  of  the  projectile  from  the 
delivery  aircraft  is  desired,  a slight  negative  angle  of  attack 
should  be  used  instead  of  a moment. 
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Appendix  A 

Typical  Wind  Tunnel  Data 


The  aerodynamic  coefficients  for  a specific  projectile  model  are 
normally  obtained  from  wind  tunnel  tabulated  data.  Table  A- I shows  the 
spin-up  data  of  a projectile  model  at  a specific  Mach  and  angle  of  at- 
tack. The  data  measurements  were  made  as  the  projectile  was  progres- 
sively accelerated  from  a zero  roll  rate  to  the  steady  state  spin  rate 
(PM).  Table  A -II  shows  the  spin-down  data  taken  as  the  projectile  was 


progressively  decelerated  from  a high  spin  rate  down  to  P > Both  of 
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these  figures  were  obtained  from  reference  11,  part  31 ■ 
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Appendix  B 


Computer  Program  Data  Inputs 


Coefficient  Arrays 

A single  000005  card  precedes  all  of  the  coefficient  matrix  input 
decks.  This  card  designates  the  total  number  of  aerodynamic  coefficient 
arrays  used  in  the  program  and  must  match  the  actual  number  of  arrays 
that  have  been  input.  The  last  digit  of  the  array  designator  integer 
must  be  in  column  12.  For  example,  if  ten  coefficients  are  used,  the 
lead  card  is* 


Each  individual  coefficient  array  is  then  preceded  by  a card  that  speci- 
fies the  desired  coefficient  symbol,  the  location  of  the  array  within  the 
input  deck  in  accordance  with  Table  B-I,  and  the  number  of  angle  of  at- 
tack columns  and  Mach  number  rows  that  will  make  up  the  array.  A 20x20 
matrix  is  the  maximum  size  of  the  input  array. 

The  first  four  spaces  on  the  card  shown  in  Fig.  B-lb  are  available 
for  the  desired  coefficient  symbol.  The  last  digit  of  the  remaining 
three  integers  must  lie  on  space  7i  9»  and  11. 
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Designated  symbol  for  the  coefficient 
r — Required  location  of  the  coefficient  array 
\ r — Number  of  columns  in  the  array 
\\t Number  of  rows  in  tne  array 


/ 11  ■'» ■ ■■ « u . , « i~» J.t .inf,, , 

! ^WiiiuniiKiJinrT*i5iiirdEnj?i^»pi^ReiEi0 

| 1. » 1 ii ..  r ,j, ,,,,,,,,,,  1 15 1 ,,, ,, , ,777..  TV  7.7777^ 


Pig.  B-lb.  Computer  Input  Card. 


The  specific  values  of  the  input  array  are  entered  by  rows. 
For  example,  the  nine  angle  of  attack  elements  that  make  up  the  heading 
row  of  the  C array  (Appendix  G)  is  entered  as: 


50.0 

" /'l«  f1  * * 'i  M **  n H'ln.iii!  «».t  » 1 rjt » « * j 


i-  « a •.  11  w mukv'Iv*  1 «i  • 


• 1 ■»  :•  1 i|<n»»iu;» 


r~o^  z.a  ‘ 57a  ibrb  15.0  zo.o ILiL— 

, , , ; j .;r» 

3.H23rri*iiiiLTrLf:D!Z!tT;:n^ 

.........  . , 

1 11  m 111  imi  1 11 1 1 > 1 i;t  1 1 1 1 1 1 1 1 1.1 1 i|n  n 1 iji  "I'  1 '"Y  " 1 ' 1 'J1 " " ' 1 1 "I" 1 1 " ' 1 1 'l 


: Ij  1 mi  mi  n' 


Fig.  B-lc.  Computer  Input  Card 


and,  in  successive  rows,  values  of  the  array  elements  are  entered}  such 


as  the  zero  angle  of  attack  row: 


*Fig.  B-ld.  Computer  Input  Card 


The  specific  values  of  each  element  may  be  punched  anywhere  within 
successive  ten-block  segments  of  the  card. 

The  remaining  values  of  each  element  in  the  coefficient  array  are 
similarly  entered  until  the  matrix  is  complete  (See  Appendix  G). 


Each  individual  aerodynamic  coefficient  array  used  in  the  program 
is  input  in  a similar  manner.  The  sequence  of  these  completed  arrays 
must  be  in  the  order  indicated  in  Table  B-I.  Only  the  desired  aerodyna- 


mic coefficients  need  be  input?  zero  entries  are  not  required.  The  pro- 


gram will  accept  a maximum  of  thirty  coefficient  arrays, 


Initial  Conditions 


The  initial  condition  input  deck  is  preceded  with  the  desired  ti- 
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The  next  card  (Fig  B-lf)  identifies  which  of  the  option  parameters 


available  in  the  program  are  to  be  used.  These  parameters  are  (Ref  5>23)< 


NORCs  Gravity  designator 


N0RC=1  denotes  constant  gravity  used. 

N0RC=0  denotes  variable  gravity  will  be  input. 


IPRNi  Output  time  increment  control 


IPRN  x At  = output  increment  that  will  be  printed.  At  is 
the  A(20)  value  (Table  B-2II). 


NALL  =1  for  this  program. 


IPUNi  Punch  option  for  angle  of  attack  (©<)  and  sideslip  (^) 
data 


IPUN=1  implements  the  option. 
IFUN=0  option  is  not  used. 


NBODYs  Reference  frame  output  designator 


NB0DY=1  output  prints  and/or  punches  cx  & ^ in  body  axes. 
NBODYsO  output  prints  and/or  punches  c*.  & Q aero ballistic 


ISCALEi  Coefficient  scaling  option. 


ISCALE=1  implements  option 
ISCALB=0  option  is  not  used 


ITRSTi  Thrust  option 


1TRST=1  Implements  thrust  equation  into  program  (See  Ref  5* 

11.  12) 

ITRSTWO  for  unpowered  projectile 


The  last  digit  of  these  integers  must  fall  on  every  tenth  space. 
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The  units  of  the  Initial  conditions  arei 


tO*  *£ 

Sec 

V V V d’  ymin 

Ft 

V V ro 

Rad/Sec 

e,  0,  V' 

Degrees 

U . V , w 

O O 0 

Ft/Sec 

I , I 

Lb-Ft-Sec 

x'  y 
g 

Ft/Sec2 

m 

Slugs 

The  initial  values  of  0 and  o< (which  is  input  implicitly)  as  well  as 
u & w are  determined  according  to  the  following  scheme: 


Pig.  B-2,  Plight  Path/Perturbation  Angle  Input 


V = velocity  vector  of  the  center  of  gravity  cf  the  projectile. 

0 b total  angle  between  the  body  fixed  and  inertial  frame. 

"Y  - flight  path  angle  = angle  between  the  inertial  frame  and  velo- 
city vector  V. 

o<0*  initial  pertubation  angle  from  the  initial  flight  path  angle. 

WD  = V sino<#  = projection  of  the  velocity  vector  on  the  body  Z axis. 

u = V cos©<0=  projection  of  the  velocity  vector  on  the  body  X axis. 
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For  a projectile  in  a 30°  initial  dive  angle  at  800  ft/sec  sub- 
jected to  a 5°  Initial  perturbation  upon  release,  with  no  out-of-plane 
components,  the  following  values  would  be  inputs 

e = -35.0,  0 = 90.0,  0.0 

u = 800  cos  10°  = 787.846 

o 

vo  = 0.0 

w = -800  sin  10°  = -138.918 
o 

Notes  0 a 90°  to  allign  the  inertial  frame  with  the  body  fixed 
frame  (See  Fig  l). 


Scale  Factors 

The  remaining  four  cards  are  for  scale  factors  (applicable  only  if 
ISCALE  was  previously  designated  as  1).  The  scale  factor  changes  the 
corresponding  coefficient  of  Table  B-I  by  the  multiple  of  the  scaling 
factor. 

The  only  scaling  factor  used  in  this  study  was  Cjj  (£w).  Cj  (6w) 
is  the  coefficient  for  roll  due  to  fin  wedge  angle.  To  scale  the  15 
degree  wedge  angle  values  of  the  coefficient  array  to,  say,  10  degrees; 
the  value  .6666  is  the  appropriate  scaling  factor  to  enter  in  the  ninth 
block  of  the  scaling  cards  (Fig  B-3). 

The  integer  1.0  is  the  required  entry  when  no  scaling  of  the  corre- 
sponding coefficient  array  is  desired. 
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Fig.  B-3.  Scale  Factor  Inputs  C(l)  thru  C(30) 
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TABLE  B-I 

Aerodynamic  Coefficient  Array 
(Adapted  from  Ref  5) 

C(I,],K)  Aerodynamic  coefficient  array 

where:  I Mach  number  array 
J.  Angle  of  attack  array 
K Location  of  aerodynamic  coefficient 


C(I,J,1) 

Cx 

C(I,J,2) 

V 

C(I,J,3) 

cy7 

C(I,J,4) 

Cys 

C(I,J,5) 

Cyp 

C(I,J,6) 

CZO 

C(I,J,7) 

Cz7 

C(I,J,8) 

Cz8 

C(I,J,9) 

C^Sw) 

C(I.J,10) 

Ci7 

C(I,J,11) 

% 

C(I,J,12) 

% 

C(I,J,  13) 

cm0 

C(I,  J,  14) 

Oily 

C(I,J,  15) 

cmg 

C(I,J,16) 

Cno 

■ecu,  17) 

C(I,J,  18) 

C|18 

C(I,J,19) 

V 

caj,2o> 

' oa.  j.  21) 

C„, 

C(I,  J,  22) 

Cy€ 

C(t,J,23) 

c* 

C(I,J,24) 

C(IJ,25) 

Olfarj. 

C(I,J,26) 

C#, 

a2 

C(I»J,27) 

cm9o 

C(I,J,2S) 

cmq7 

C(I,J,29) 

Cmqg 

ca,j,3o) 

% 
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TABLE  B-II 

Computer  Input  Array 
(Adapted  from  Ref  5) 


A(l) 

lo 

A(12) 

”0 

A(2) 

xo 

A(13) 

vo 

1 

A(3) 

yQ 

A(14) 

w© 

A<4) 

zo 

A(15) 

H A A 

| \ft# 

A(5) 

Po 

A(16) 

V 

A(6) 

V 

A(17) 

g 

A(7) 

ro 

A(18) 

m 

A(8) 

0 

A(19) 

d 

I 

A(9) 

<f> 

A(20) 

At 

9 

A(10) 

* 

A (21) 

*E 

A(ll) 

1.0 

A (22) 

66 


AFIT/GA/AA-77D-8 


Appendix  C 

Computer  Output  Example 


The  computer  output  parameters  are  shown  In  Table  C-I. 

TIME  is  measured  from  release. 

RANGE,  ALT,  Z correspond  to  x,  y,  z of  the  inertial  reference 
frame. 

V s total  velocity  of  the  projectile  with  respect  to  the  air 
mass. 

p s spin  rate. 

ALPHX  .,/«*+$*  (where  - -r  denotes  vector). 

M s mach  number. 

PHI  * roll  angle. 

ALPHA  = angle  of  attack. 

■ 

BETA  * yaw  angle. 

The  next  set  of  parameters  are  useful  for  stability  analysis  appli- 
cations.  The  outputs  are  based  on  linearized  equations  for  the  tricyclic 
response  of  a projectile  to  an  initial  disturbance  (Ref  12,  8 or  11 )i 
L-N  e An=  damping  factor  for  nutation  mode. 

L-P  = Ap  = damping  factor  for  precession  mode. 

W-N  =CJn~  nutation  frequency. 

V-P  =(Jp  b precession  frequency. 

S b s > stability  factor. 

TAU  b T ■ dynamic  weight  factor. 

K-T  * * trim  arm. 


AFIT/GA/ AA-77D-8 


/Q  PPEN  Df  X D 

ROTATIONAL  PARAMETER  TRANSEORMAT/OAJ  ( REF  10) 


zy.  y.  s INERTIAL  FRAME 


x f K « INERTIAL  UNIT  VECTORS 


xvz  = body  nxe’D: 


TUE  UNIT  VECTOR  V /S  MADE  UP  OF  THE  THREE 
COMPONENTS  Y.  ft  ft  50  THAT 

V = V,  * +Vxj  + V 3 k. 

NOW  DEFINE  V.  - X,  / 5/A/  ®V:z 
ft*  W 5/ A/  c^/2. 

Vs*  Xs/  SlN<*/2 

where  x.  Xi  Xi  /?/yz?  <=k  /?/er  UNKNOWN  parameters. 

USmCrTUE  REIjOTION  FOR  DIRECTION  COSINES  TV/FT 
THE  SUM  OF  THE  SQUARES  OF  THE  COMPONENTS 
OF  A UNIT  VECTOR  MUST  EQUAL  ONE  C-IVES 

AT  +-  x\  - X^  = SJN^o</CL 


1\)d 

LDo 


£3)p 


2.  2. 

A/0U/  DEFINE  ANOTHER  PARAMETER  \o  - COS  <*/2-  , 
SO  THE  EXPRESSION  CAN  BE  WRITTEN  AS 


X.  *-  2.  2- 

A/  + Xlt  + X3  + Xo  “ I 


( 4/p 
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THE  BODY  FIXED  FRAME  CAN  EE  RESOLVED  INTO 

inertial  components  as  follows: 


X1 


* 


LET  /}  BE  A UNIT  VECTOR.  ALONG-  OY  AND 
& ALONG  OIL  SO  THAT 

A - CK,  .T  a.  CL*.  f + CL  9 < ' 

B = b,  J -f-  bx  £ 

where  a . ci i a~  , b,  b*  bn  are  i/jeetia l frame 
COMPONENTS.  ALSO  LET  ^ BE  THE  AA'.QUfJT 
OF  ROTATION  OF  THE  PROJECTILE  ABOUT  THE  OX 
AXIS  WME/ZE  r IS  A VECTOR  FIXED  ON  THE  BOOT 
BEFORE  ROTATION  AND  F'  IS  IN'S  POSITION  AF1ER. 
ROTATION.  THE  DIRECTION  COSINE  MATRIX  FOR 
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7 11/ S Z.OTAT\OU  I £ 


I o o n 

0 C&t  rSe*  n 

o s-i  r*a 


U*)p 


- Fx  £«*  - Ta  S<=-< 

, (7 

f/  - r-x  s<=*-  + r3  coc 

V/HEPE  r>  Ft-  Ei  ARE  COMPONENTS  OF  t~j 
C = COS  , 5 = SINj 

AJFCAT/VE  POTATO  AJ  ABOUT  X U/AI  UiO? 

573  EON  1 2-  SIMPLIFIES. 

EXPRESSING  fs  COMPONENTS  IN  TUFT  INERTIAL  FRAME  . 


A a.  ^ 


A* 


^ Y “ siN^/z  ^ + 5 W^/T  ^ + SIN*/  2. K 
ka  - a,  a -t-  a f & 3 & 
rsB  = b.  a + b -t  y 3 * 


COMBINING-  EON  7 * S’  IN  TERR'S  OF  THE  BODY  TRAWL  l 


'7  A.  f * /u  * -+. ./.1£. 

»•  » * 5 IN’YiJ'  SIN*-/!  1 SlN"</a 


Al  i*  -4-  £ 


r*4t(n,  uv* 


ui  $** 


A 

- b.  £»*-)  .<  -*■  ( «a.  5'*°*-  * *>*•  c0-50')  •/  * ( * a -'■ b3  <ro-2^  K 


1 


DEMOTING  WE  TRANSFORMATION  MATRIX  BY  [c]/  THE 

exrhessioju  Tim  trawfoltis  we  components  of  the 

MOVING  CODY  AYES  SYSTEM  TO  IAJERTIAL  REFERENCE  IS: 

“1  r . *1  r — 

C„  C.-l  Cfz  ~~r,^  a.  b>  * a.w-4,c« 

c..c^c„  a,  b»  - -£^,0,01- I»s«,  o.s^^fc*OK 

Cz.  C,z  C»  5^  a>  bj 

— — — **  *—  “ 

. ca/v  (/o)p 


POST MUIVFEDNC,  EV  7WE  INVERSE  ( TRANSPOSE  FOE  HIE 

okti  io  qoajal  syste:  /Vi ; : 


w.  Sin  o'/a.  ' 

£ii  Cm.  Cli  — 1 

C..  Gr  G-i  5^^ 


— 

V 

V* 

SIH^h. 

Sin  »Vj 1 

SW'A’z 

a, 

ax 

ay 

b> 

b*. 

bj 

EQM[/0l 


MLT1PLY/N6-  IT  OUT  I 


Cu  C.t.  c.  3 £T5£  + f «N  ft ) c*  , + (a.  b -dib,)&  + (*.  <?*  + b,  b» ) c-< , 

7^4  fa.  i>3  - b.  <3j^  S**  ■+  ( b . bj  -+  «.  ic°^ 

G'  Caz  Caa  “ -h(b,aj  -a,k.)S*<  y 'pi/E*  b’JC'x.  > 

+ -+  (azai+  b,ba  ) <?<* 

« 

Cji  Cm  Ga  hi}ll  + (tt,a3  + b,b3)Gx  +(biai-0ib3)  S+i 

s2-'/!  V-2- 

j~_^+  lojdi*  b j bi)C<*  +(l>i03-  bidi  'F°t * 


£fOr/(/2)p 
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THE  FOLLOVJINO  ZCLAT10/JSUIPS  UE7 UE.E/J  DIZL'CTIOJU  COSINES 
APPLi  TO  AN  OIZTNOCQNALJ  Zl&U'f-UANQ  sry STEM  l 


at"  ■+  a\  •+  a\  — I 

EH 

.2-  .x  .x 

I • bi  •+•  bx  ■+  b 3 - I 

a,  x ».  ■+  $ 3 x ■a  — o 

b. X.  + b^xt  +■  bsXa  = o 

|<3i  S,  -*  a-L  bz  + flj  bj  =.  o 

b.  = — - Ajd.) 
b.=  sW,(A^,  - x,  a3) 
bo  = L/j ( A.Q-l  — Ai  £() 


X, 

X, 

X 3 

S-VJL 

S~*/x 

a*. 

a3 

b. 

b*. 

bj» 

(n)0 

(M)d 
<JS)o 
(.11.)  0 
0 7)  o 

0s)o 

09)d 
(20)  0 

(2i;c 


I 


i i 


EQUATION S>  US)  AND  UU  ) CAN  C£"  h/PJTTlzfi  AS'. 

CX-x,  ~ - ( Cl,  >,  + G3  Xj)  (22) d 

ba  = — ~ ( bi  A , -t-  hi  \i)  (23) o 

SUDSTITUTINCr  (ZL)k  (22)  INTO  L 17) 

cr.bi  ■*  tfi  hi  -t-  (3j  hi)(  b 1 hi  bj  AjV  *7.2  £>j  = 0 

(2-1  )d 

a,b,(  xl,  * X*)  * ^3bi  C A\.  + &)  + (03  U •»  x,  \3  =<2 


I 


i 
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SUBSmUTINt  (22.)  INTO  (IJ)  AND  (23)  /WO  (/-?) 

0*(  A*  ->  >x»)  r Xx, -i-fj)  - X.  - 2 A. A>  <7,  <7  3 (25r)p 

bif  A« ■/  Xi)  >+  ba  f Xj  + Xj)  * Xv»  -2  X# Aa  bi  bj  (2£>)p 

g-QU/IT/ONS  (22)^23)  INTO  (2  I ) 

-hib*  + \zb,a3  - *2^7 

4 tfjXj)~Xa  a,  + hg(b,)u  G)* 

/)F7££  CANCELLATION  S AND  gEC-fZOUPlHO 
-bj£T,  [ X%  >fi'  XsJ  ■*■  b.  a 3 [ X/  ■*■  Ai  * X3  ] - XjS'T'-a. 

MO  SINCE  X*+£l-»  X?*3  - I 

b3'-S7  [b'<*3  - Vs~/J  (27-)o 

EQUA7/0A/  (27)  //VEO  (2^) 

a, b, (&*&)+  |j(X^X%Xb4a-  - A.Xaf2tf  »b.  - - 0 

^?7b.  {£+X*J  + * °j[1XiXj  tf'b*  r ^^A.AiXj 

EQVA7J0AJ  (2.7)  JWrO  (2t) 

£ (*,-£>-£  I bA  -Jy’uV  >>M  = x\  - 2 |r(  b aj!  - r-U  *•  x 3 

a!b'.  (M.)-dj I>.MVX!j)  + flab.  [2-Ji.Xj  o.f.  - jb^i  if**-*)] = 

A*  d?  + lb, a,  - -^(Ai  -Aj)  + *»b.  ^MX\^) 

LEFT  S/DE  OF  (2S)*  b,  - LEFT  S12E  OF (29)  TAElLEFOfZC. 


U«)p 


(29), 
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.|^A,AiA3  = \\aT-+2b,a, 

mLVPLV/Nt  n-m  cv  -^i 
da  - fx\ *x'j ) - £ -ya  Xx^jr,  “ tf.X.Xa}  *r; 


,.*Aj 


(30  )p 


EQUATION  08)  AfiO&'f) 

b,  ♦ (&*  £-J*> 

eq  uati an  (?o)  into  (33) 


(3/)d 


b,  = + 


s**  x? 


(Al.  + A-  )[-J. ...  )f[ A*  (£*>'>)-  tfxi  s~/J±  - a,h 
U*yi  5^x4  x^yjLs-'/i  Jt*  J 


-- '&£&  - r'/JU  - 


EQUATION  (32)  IS  US'ETD  TO  EVALUATE  C„  IN  (/ 2) 


(32)o 


C"  = + ^»v+^>C6s<x 

FROM  THE  DEPIIUI1IOUS  C05^/x  AND  COS 3*<-  2.COS*°c  - 1 
Cn  - + {Xv+Xa)(2X<»  - l)"j 

! — fx2;  ♦ 2 X*  - / - 2 Xo  - X»  - 2 A1/  X%  * X,1 

• "5v-/xl  -j 

• r _L_[2a%Al)  - 2fXv.*-X'lXlo  - Xo-l] 


?5 


= 2 { X.  ■*  >0  )(l~  \)>  ) - ( I - }ia)J 

= (LrJp  foot*  \l)  - l"] 


usmcr  s1-'/!*  cWx  s /-  x; 

ft 

C„  ~ - 1 

(X3jD 

• 

sz/mil^rlY 

1 

Cxx=  4-  (#1+  b \)cos*t.  = 2(x\+  x*.)  - z 

— — 

(3^)d 

■ 

1 

£33  - £ )cos~*-  » 2 CA3}  -X*)  - | 

(T5)o 

) 

1 

FROM  EQUATION  (2.7) 

Ml 

b.tfa  - b3  = -Ai 

(^)o 

BY  ANALOGY  AND  UITH  ( 2-0 

■ ; , 

birf,  - b.  a,  --  ^ 

(37)p 

!; 

ba^x.  ba  i * 5 ot/i 

(3S)0 

evaluate  a,ch  + b,  b,  by  using  (27) 

+ b.  bi  - <7i  <7i  + 6.  -i-fbi 

= U'^-ara 

EUB  EQUATION  ( 324  *. 

- ^3  Aa^Xj  _ bi  X 1 

~ a,  s*-*Vx  a,  s-vx 

SUB  EQUAlld/J  £30): 
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: s^vTa, f ^ A',J ~ S-KK,  a.’J-a.x.Xj}- 

— Xq  / v ■»■  ^ v1  \ __  X?  fit  _ X»  X 3 _ bi 

” SJ*/i  4,0,  x J r-Vi  b,o,  rr*Va.  a,  bi  S-Za. 

53-Vza,fa»  S'Vi^ibr  S3 


SUB  EQUXTIOA)  (32.) 

a, ax  -t  b,  b2  ■=  - -A;-A?- 

.Sz  -Va 


SIMILARLY 
bi  bi  — — 

a2Oi+  bx  b3  = — 


Xi  Ai. 

Sr,c*/i 

X r X 3 

S^OL/Z 


EVALUAT IM6-  Cn  FROM  ECJnOI)  WITH  (-ftf)  »/YP  (37) 

£Tn  = jx^r  +(&i  Gi . ■*•  bi  -f-  bx  — b/)  SJA/oL 

- -AA-  _ _2lj£x  (2  >£-,•)  4 r-/2  xe 

■S^Wz  SMi  ' $~/z 


WHBF-E  COS  2(-7i)  - 2 ( Co^-vi.)  - | 

5 W 2 (-vi)-  2 ( SW*/*  X Cos  -"/a) 

£<3.  = 2.[  Xi  Xi  -t  X0X3] 


ccs“yi  = Xo 


(39), 


F°)t 

Me 


IN  THE  SAtVE  A/IAA/A7ER.,  US  IMG-  EOUATIOAJS 
(32)  (37)  (3  8)  (3  9)  (-7-0)  (-7 /)/  77-7 /T 


A FT  T/ OA/AA -77D- 8 


REMAINING-  COIAPOfJENTS  ARE 


Cx,  = z[ 

A»  Xa  - An*  Aj  ] 

(-«)0 

C,  3 r 2.  [ 

A#  A3  — A<j  As.  ] 

c-w^ 

Cn  - 2[ 

A 1 A3  — A a A 2] 

Cn  = Z[ 

Ai  Aa  - Aa  At] 

('Jfc)p 

Cjz  = 2.  [ 

"A*.  As j — A#  At] 

W„ 

SO  Tit  AT  THE  DUIECTTJ/J  (ZOSIAJE  MATHl'K 
BETMEEKf  IV E INERTIAL.  AND  .BODY  PAX  ED. 
REFERENCE  FRAMES  AFTER  ROTATIOAJ  <=*  IS 


Co  Cl  X ClT> 

2(X^Ai)-  1 

2(  A.Az  + AoAs) 

2/  A»  As  - Ao  Aa  j" 

Cxi  Cxx  Ciy 

— 

2 ( AiAi-AoAa  ) 

2 ^ ??:.+  A>)  ~ / 

Z( AzAj  + Ao  At) 

Cji  Gt  C31 

_2(A,A2+AoXz) 

2(  AzA.s-AoA.) 

2.  f A 3 f-  A\ ) - 1 

— 

E(3N  »)D 


WHICH  IS  FQLnVALEAJT  TO  THE  EULER  ANGLE  FORM 
OF  THE  DIRECTION  CO  5 1 ME  MATRIX 


oo  cv  se 

-SV' s<p  - CWS9  c<p  CO  C<H 

-SVCQ-+  CTS9  2<J>  -C9S<t> 


SW  C& 

CV/Trp  -SV'50c.(p 
CVC6  + S¥'S£/  5<fl 


EQN  09)d 


i 


?8 
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l 


CAMPAEIN6-  COEFFICIENTS 


£//  = 

if  aUaj 

- / ' = CV  C9 

(SO)o 

£ix  - 

2 ( A*  + A\} 

- 1 - C <t>  C6 

Uri)  P 

C33  = 

2-(  A*  + A^) 

-1  - C Y>  C<2>  + SY'S&Sd) 

(-52?  o 

FROM  EON  U) 

. 

& ^ 

1 - c ?u  + \h 

- + A*) 

to)p 

EQLMT/ON {S3  ) INTO  (S£>X£-lXS2.) 

1 - 

ELi  A>  Aj)  = 

CW  CO 

U^?d 

/ - 

11 

8 

+ 

H/< 

rJ 

Cd)  CO 

(52r)p 

1 - 

2 C A* -*  A*J  * 

CVCQ+  ST'SOEcp 

c**i?p 

EGLN(£T)-  (£5-) 

* cvccd  +srsos-([>  -c&co 


EON  UK)  + (57? 


l-4^=  CV(CO  + C<t>)  + EV'SUSV' - C<D  CO 
4 A*  = / -C.V(.Cd  + C<D)  +•  C(0  CD  -SVS9S<p 
\J£INC  TI11C  IDENTITY  C&+C<D  ^2  C(jt±P)  C(£=*) 

I - 2C7/[  C C 2^]  -t-  £6)  Cd)  - SVSQSCj) 
USIND  THE  IDE/'JTITY  C 2.r=cV  -sV 
4 Aj=  I-2(CZ  - 5^)(K|  c*i  _ Sj?  ]+  CQC9-  WSO£(t> 


4\x-\+2 

-2 


1 2.  X 

5 ® 

1 2.-.  2 

£*V  £ltf  /*■© 
Z I 2. 


1 111 
+ SZCM  C V] 
2 i 2 J 

-+-5V  s’*  slsl 

2-  2.  zj 


C<6  C£?  -ST  26 50  + 


vsinc-  identies  a 2.0-  zbo  - / 


«£  C26)=l-2sb 


(S7)  o 

(SS)0 
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4\  - 2 fc's  S*J|  S£.  -+-  Si£  C£  C\]  -*■  J + 

- *[(C v+iXcoiiXc4+0+0-cvXt-c9Xt-cd%  + 
+ C(bCQ  - SVSQSCS 


s\  -tSZC^C^]  + /-«- 

- ^orcd + c vc  <d  •*■  cdcv  + /]  + 
4-a$ce  - sv'z&sd) 


"--if 


<f<£  S'S.SZ  4-  5V  Cjg  cVj 


L 2 


-i-j^|  — CV(CQ-*-  C-(D ) I-  CtVCQ  - SVS6)S&] 


_ ’svzescD 

z 


f 

4X\oT&  S£  S*®  4-  sVcif  cV|  +2M  - 1 SVS9S0 

’ L.  2-  2-  2-  2-2.24  2. 

U5/N6-  TW£T  IDENTITY  S D.Q  ~ TlSV  C& 

SB. SSL  4-  S J'  C*g  <f<2  -2  S£  C|  5|d|  S-| 


2 2 2-  2 2 
V/HICH  F^CRJZS  INTO 


)\x-(CXSRS&  - SZ  C§.  CffiV 

' X Z 2 2 2 2 ' 


7S  r C vi  5££  S£  -5J£ 

2 2.  2 2 2X2 


B7  ANALOGY 


- - £(0  S £2  S5£  4-  S &C.Q  C<£ 
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EQN  (n)  -(S'?) 

4^-1  = cvcb  STS&£cf>  -C<t>C-Q  - cvco 
4 / —C&CCV'+C<t>)+  C&  CV  1-  SVS0  Sep 

rue  souu  to  Tins  \zqn  is  similijwl  to  izqn  s ~8  : 

■ I s*  * cf  Cfci) 


eqn  (s'*?)  y?vro  (s~3) 


i -T £.*±(1-  crc0))  - -if /+  cv'c*)-  x* 


4-2C&  S & C.SS  SSSCH.  s& 

Z Z Z Z Z 2. 


fetfi+crce)- (i-sfid&dq  -s\dq  dq  + 

+zc-&.<±&<Lr  siiL<L&.<a 

Z Z Z Z Z 2. 

Aio-i[l+(2C^-,X^CZ^-/j]  -rdqdqfz  4- 

-£JZ  <,V  _ s\6  £*t?  +2C5 5^  CV'Sc^  C9.SE 

2 2.  2 2 2 Z Z '2L  Z Z 2. 

A\=  I +■  C\£  (I  - ) -+-  /s  sk  4-  c*V  eda  + 

° 5_  Z ' z'  Z Z i 2.2 

- c7^  - +.  ZLC&.S&  C-SSS'JSCOSO 

2.  2.  Z 2 2.2.2.  2 2-2. 


nots  Twsr 


i + css  c’a  - css  - c%?  -do.  dss  - 

Z Z Z Z Z 2 


I + dQ  t - dsL  ) - cV  - 5^  skf  = 

Z Z 2.  Z 2 

/ - <rk  - cdq.  ~ o 

THE(ll£F0(l£ 


•>  XX  X ,. 

X,  - £ t?  C.<£ 


-t-sk  s\g  sk  -ricdisdi  csssv'  case 

Z z a.  z i-  z -a_  z 2- 


A,-  CJH  C&  CQ.  -t-  SXSfaSR 
z .1  z z z Z 


M) 
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Appendix  E 
Coriolis 


From  the  klnematlcal  expression  for  motion  in  terms  of  moving  ref- 
erence frames  (Ref  7*111)* 


(1)E 


where  the  superscripts  denote  the  applicable  frame  of  reference. 

The  term  2 CJ*1  x r is  the  coriolls  acceleration  and,  for  the  ap- 
plication of  a projectile  moving  with  a velocity  (V)  over  the  earth  ro- 
tating at  rate  (/l),  the  coriolis  expression  becomes 

2 A x V (2)E 

Define  the  body  frame  (XY2)  moving  globally  with  respect  to  a north- 
east-down reference  frame  (NED)i 


Fig.  E-l . Local  Level  Reference  Frame 


where  is  the  heading  with  respect  to  north  and  L denotes  latitude. 
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The  transformation  between  the  body  frame  and  north  frame  is 


COSY  SIN  Y O 

- SIN  V COSY  0 

0 O I 


(3)E 


The  earth  rate  component  as  seen  by  the  body  frame  is  then 


COS  Y SIN  Y 0 

~ SL  COS  L 

= ci  rt  -- 

- SIN  Y COSY  0 

0 

0 0 1 

-XI  SIN  L 
_ — 

neon  cosy 

- ILCOSL  SIN  Y 

- XL  S' IN  L 


(4)E 


so  that  the  coriolis  expression  becomes 


2.-0.*  V = 


2 


yt 

SL  CL  CY 


1 

- XL  Cl SY 
Mr 


K 


-X1SL 

V0 


(5)E 


where  V^,  Vg,  VD  are  the  N,  E,  D velocity  components.  The  scalar  com- 
ponents of  coriolis  acceleration  then,  are 


(2j7xv;y 

(2J 


-2VpXl  CLSY 
- 2 Vp  XI  CL  CY 
2\(J1  CL  CY 


+ 2\(XL  SL 
- XVkXLSL. 

+ 2^X1  CLSY 


(6)E 
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AFIT/GA/AA-77D-8 


Appendix  F 


srxoor  TK/Tk  opt«i  hh  «.$•<.!■>  oi/i t/rr  u. 09.00 


t PROGRAM  SIX30F  ( I MPJT,  OUTPUT,  T»«*E5*I  NPUT  ) 

EXTERNAL  APOCNA 
CALL  CFRO 
CALL  I-IIT6 

s cut  sixoec  uRtnmw 

iro* 

CN3 


AFIT/GA/AA-77D-8 


limOU’tME  A. OCHA  74/74  OPT  * 1 PTH  4.5»414  IV1  7/77  17.09, 


1 

SU9*'HITIKc  AROCHMH) 

C H0B9.E  ?'C  9.  S.  FtK‘-M17  <°r,  «**?0-SPA2F  FNG9. 

W9LC12G6 

c 

UNlVFRSITr  or  MOTPc  Oil?,  *OTRF  9AMF,  INO. 

M*l  0 1 ?r*7 

c 

ML0t?5e 

s 

C IPO 

C 195*  MOPFL  ATMOJPMFRS  SU9R0UTfNE 

M9L01259 

c 

30HMOM  /iIR9L«C/  T#®,o,V<,»/5 

W9L01260 

M9L01261 

oiPr,iricN  a<3i  ,hji, ci»>,iii> 

W*t  01262 

01 9f 931  ON  H9(12),H1(  1 1)  ,47(11)  ,'43(ll>  ,TP(ll)  ,P9(M>  ,P9IU> 

M9L012S* 

It 

c 

W3L01264 

3AT1  M9/C.  ,110  00.  ,750 31.,  47  0 11. ,53000  .,/  >00. ,90  000. , 105  000.,  160  00 

00. ,170100. ,701003. ,7030  03./ 

DATA  Kl/-.  7755617--.  , C . , . 1 ’ 94t67-« , 3 . , - . 1 53  70  2F-4 , 0 .,  . 751 55  IE-6  , 
t.05S759E-4,  .75.1.  1--5,  .’5071  tr-5,  .7771  ?i-:-5/ 

IS 

DATA  M*/-.5?56l?rl,C. * .11 7 91’r?*.C,-.759?13El,0., .55417071, 
I.1T0174F1,  . 341 641* 1,  .fti’Tli'l  ..  3761  3701/ 

34TA  W3/0.,.157S<9r-?t  7. , . 12*5fc*E-*,  C • • • 23 6? *4P-3 « 0. « 0 . • 0 • • 0 • • 0./ 
DATA  Tl/. 51069  05  3,.  7*31.9-)-  J,  . 1199117 3 , .50 1 7 95F3, .5 0575 "73, 

t. 7951 51 F3,. 79* 151 7 7,  . 4 15!  ■ *r  7,  . 2 *361  5 74,  . 776615E4,  .75361951./ 

70 

DATA  p9/. 211671774, .4T'‘7*73,.5137Q7?,.?5155El,.12151El*.2t  057-1, 
0.210  397-2, .155677-1, .75  7«r-3, .50 9547  - 5 ,. ’97597 -5 / 

DATA  93/.73T6977-7, .70677-1, . 7’657- 4, . ’312 4E -5,. 1394657-5, 
0.411*97-7,  .47617-5,.  ?7  7?r-9,  .114  57-1  1 , .1  339E-U,  .61137-17/ 

DATA  CONI , CON?  ,C1N3, 30  *(4,0095, 2096/. 3041,  5356766. ,49. 070576,0. 0226  W3L01771 

75 

t355F-05, 191.77,9:033./ 

M9L01279 

0»T«  A/1.,. 751511, .3*5717/ 

OATt  9/0. ,.174164, .773966/ 

OATt  C/0., .77E6,. 1176/ 

DATA  0/0., .25E5, • 14F6/ 

30 

c 

HGP»20Ni*M/Cl.M30Ml  •H/COM2)  ) 

M9L  01202 
W3L01283 

IMHGP.LT. 0.)  MG®»3. 

M3LC1294 

00  1007  9=1,11 

M»L01?*5 

1MH‘.P-H?(M)  ) lQOTflOC^tlCO? 

MIL  01256 

35 

1007 

DONTIN'JE 

M9L01237 

IT  <<MGP-H0C12> > .GT.Q  .)  GO  TO  1052 

M0LO12* 5 

H*l? 

W9L012S9 

too* 

H»M-1 

M9L01290 

1004 

TMsTI(S)*  (t,*HlM)  »CHGP-H3(H>  )) 

M3LC1291 

>0 

I r C( ”G P-90300.). GT .C«)  GO  TO  1006 

MHO  129  2 

T*TM 

MHO  129  3 

GO  TO  1070 

M9L01294 

1006 

IF  ((M'-,P-1S9COO.>.GT.O.>  GO  TO  1009 

MU  01 29  5 

1*7 

MH01216 

05 

GO  TO  1007 

M9LC1297 

1009 

1*3 

MU  01 29  5 

1007 

T = T9*  (A  (I) -3  (I)*  AT  AMI  (Hr,o-C(I))/3(I))  ) 

M9L01299 

1070 

if  <h?(h))icii*i3?o. ion 

MU  01300 

50 

10  11 

TPMPs 1 • 4M1 (M )* (HGP-HO (M)  , 

MU  01  301 

P*t“J(H>/TCHP*#W2<M) 

W9LQ130  2 

RxPq  ( M)  / TEMP*  • (1.  M2<»1>) 

MU  01*0  1 

GO  TO  1030 

MU013C4 

1070 

TTH:>*rxp<-w3M>  •<m;p-ho(i)  ) ) 

M9L0130S 

55 

P=PB(*)*TEMP 

MU  01  30  6 

Ps99(M) f Tfl® 

MU  01*0’’ 

1030 

tr  CIHGP-C0M6I  .GT  .0. > GO  TO  10*2 

WH  01*0* 

V9*rOM3*SORT (TM)  • 

MU  0 1 309 

V«C*WK»*  1.5/C  CT»CON5)  *P) ) 

WU  01310 

50 

?rTU*H 

M9L01311 

1052 

r = 0. 

MU  31*1? 

PxO. 

MU0131 3 

R»0. 

MIL  Cl *14 

10  32 

vs*c. 

MU  C 1 31  5 

65 

¥<*0. 

MU  01*16 

f FTUPM 

mu  fit  *17 

*MJ 

MIL  01319 

BEST  AVAILABLE  COPY 


AFIT/GA/AA-77D-8 


f 

* 

i* 


timtrii’  c»«o  rs/rs  opt«i  pin  s.5»»ii.  ovi ?//?  ir.ot. 


i 

t 

C 

suo’ouriNE  cfro 

00003070 

00C93CA0 

C SXlDEG  SO  C.  «.  INGRAM,  afrr-$pac*  MG*. 

00000090 

e 

UNIVERSITY  OF  NITRE  OA -»£! . NDTRE  91HE,  |NO. 

00000100 

> 

c 

03090110 

V 

c reads  cor*f ieirnTS  in  tabular  fo°i  as  a function  op  ingle 

OF  00000120 

* 

C ATTACK  AND  MACH  NUNRER 

00000150 

C NAXlNUM  MUMPER  OF  An;L"S  Or  ATTACK  = ?0,  NAXliUM  NUNPFR  OF 
C NUN9ER5  * 20 

Jt 

o 

X 

tt 

C 

CQQQG15Q 

• 

toot 

FORNATC  A4.1X,  312) 

000001S0 

100? 

FORUT  f //l 4X , • At  3H  A * , ?X , 1 0 r 1 9. 4/2 1* , 19*11  .4) 

too* 

FO?N.VM/5X,*HACH*  , IF  lO.fc.’X,  1 9r 19. 4/ 21X, 19*10.4) 

10  0* 

*0»HAT ( /1XA6, 3X2HC(I2t  IN) ) 

00000190 

i 

is 

10  06 

FORMATCSF10.4) 

loor 

FO*HAt (detO. 4) 

lONHTN  /COEF/  «MCH(?r,  TO)  .tL’HMEO.IO) ,C ( ?0 . 20 , JO ) ,NO» (30 ) ■ NOM C J01  00 0 00 ?0 0 

lOMNOII  /COEP1  / N°L  Y 

09000210 

COHMtWI  /POL/  CE(YO),  »L0HE»,r.'i»,CMPA,C7» 

0000C220 

i- 

n 

COHYON  /SLK1/  R'JY 

OIHENSION  AHC(600)*ALP(SGP)»'C!1?0P9) 

0000C2T0 

EQUIVALENCE  (AHC,  AYCHt  , (ALP,  !L=-IA)  , (CC,C) 

00030240 

' 

c 

000  OC  2*3  0 

■ 

\ 

C SET 

ARRAY  TABLE  TO  I*  R9 

00003260 

» 

N®L  Y*  0 

00009270 

90  110  1=1,600 

00  0 OC  29  0 

: 

ANCIT)*0. 

000  10290 

• 

■■ 

110 

ALP«T)  = 0. 

00003300 

' 

90  111  1=1, 1 ?0  99 

00010310 

\ 

31 

111 

:cm  =0. 

00000320 

00  11?  1=1,30 

00000330 

NOAH  1 = 0 

000093**  C 

] 

NON (T 1 * 0 

000003^0 

11? 

CONTINUE 

00C003S0 

35 

CZA.P.O 

00090370 

i 

CNA.0.0 

00  C 00  33  0 

\ 

CNPA*  0,0 

C0109390 

c 

•0000400 

u 


<1 


55 


St 


C REAO  NUH9FR  OF  AEROOTNAHC  COEFFICIENTS 

000C0410 

1005 

RFADf  5,  1095)  PUN.NC 

F09NAT(A6,4XtI2> 

100  0 

PRINT  10C3,RUN,NC 

FO»NATf 1M1 f 2X, A6, 2X, 12  //  ) 

t (O’ 

CONS* ANT  ANGLE  OF  AT  TAG*  NO  ANGLE  OF  ATTACK  TABLES  ARF  READ 

00  0 09430 

C FOR 

CON  T A NT  KACH  NUN9FR  ONE  T A aL*  OF  NACH  NJM9ER  IS  READ 

00000440 

00  1*0  M=1,NC 

00003450 

REAOC 5,  10019  NAMrf  <f  »J0A<<)  , NONCO 

00000469 

PRINT  10  04,  NAME,  K 

09000470 

NN*MOMCK) 

00003430 

IF(NM.EO.O)  NM*1 

00000490 

NAr NO A ( K ) 

090C0509 

IF(NA.FO.O)  NA  = 1 

•0000510 

IF(NA.FO.l)  GO  TO  119 

RPAOCfo,  10  06)  ULPHA  C I ,P) , 1=1 ,NA) 

PRINT  l 0 C 2,  CAL°<A(I,K),I=1,NA) 

00000540 

lit 

DO  1?0  1=1, NM 

1?0 

R*Ani5,  1007)  ANCNII,K),(C(I,Jt>r)9JrlvNA) 

t?5 

90  12*  1=1, NH 

00000560 

PRINT  l 005,AHCM(t  #<)  t CC(1,J,0  f J*1,NA) 

ooo  oos^o 

no 

COMTIM’IC 

00003580 

RETURN 

00009590 

END 

BEST  AVAILABLE 


i 
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AFIT/GA/AA-77D-8 


suMmrmr  to 


7k/7k  0°T  «i 


rTN  4.5*4J4 


09/1 2/77  17.09. 


X**V*f7X,*°*,3X, 

•L-N**3X,*L-P*f 

♦4X**P*/StC*,2X, 

,*i/stc*#ixf 

F6.2, 2F7.1, 


SIJ9P0UTINE  £0 
Rr  A(.  HAS5tXX,lY,L 

XHTfr.ro  p g 00004990 

C0K-9M  /AIR91K/  TE*,PRCS,R>0,<VI>,VA  00004600 

30»W*»  /COCF/  A m:m<  ?•:,  3G1  (?C,XC>  ,C«  23 ,20.70)  ,NOA  HO)  * N9H ( 3 0 > POOC^SIO 

50MMW  /IN99LK/  I 9 c N , I WIND* I Tro,  IDRN,  I ALL  , N0RC  , X PUN,  N909Y  * I TRS  T , NT  00  C 0 -3  20 
HST  00  0 C 46  3 0 

SOW^ON  /INIT/  TrrL'(12»*SCA:F(33l  f A(?2)  00004640 

DONMVI/TMRUST/TlNCtS)  ,TpsT(1?>,NAS9(16),CG(15> ,IX(15> , IY  (15> ,CGP( I00JC4650 
Si)  00004660 

30M*nN/TST0AT/FO$,  T-l  YC  w ,*L  00004670 

/H1N3LK/  N94,YW<1?9>,W(1?3,3>  000C4630 

20H99M/C0t®l/  NP.Y  00004690 

S0MH9W  /POL/  CE  ( 3 1)  , A L°W A= , CMA , DNOA , C*A  00004730 

30**M9M/  F09K/9t(14l  ,37(14)  , I- 1(1  41  ,3'  2 (141  ,WlN(3>  ,OAL(10>  »39t(10>  00CC4710 

COM^9*l/tC9K?/Y,0,  AIRMAN,  A“C4*I,  PG.LN  00004720 

COK*n*)/E09<3/AxX,  AXY , AX7,AYX, AYV,AYZ, A7X,AZY*AZ?  000  04730 

30WN9N/E03K4/  LC3UMT 
50KN2N  /°L  Kl/  R»JN 
9ATA  00*10/57.295779/ 

1900  FO^i* T(56Xt*--r*> 

10  01  FORMAT { IX,  T4  **TH*%  5X,*P«*r,c*f6X#*4LT*f  7x,*Z*,9X,  *V* » 7X , 3X, 

1*AlP**A*  t3X,*H*  ,5<  ,*p mi  #,?*,•  tL® HA  *,?X  ,*r,"ra*,2X,*t-M*»3Xt*L-P** 
24X,*W-N*,4X,*W-t>*  , 5X,*5»,5X,  *Ta-J*  ,3Y,  ♦<-r*) 

10  0?  C07M*  T( IX, T4, *SFD*,7X, ♦cP£T*,6X,*PFtT* ,5 X , • eF£T* , 4X *• FT/ScC# f 2x , 
l*9A/5FC*,2X,*0fG* , 1 3 * , * DP  5* , ? X, • D^G* , 4 X, * 3F 3*, iX  , • i/StC* , IX, 

?*l/SrC*  ,1X,*»A/S£C*,  1X#*RA/Sr5*  ,tSX,*OFG*f > 

1003  FORMAT (ix, lPe.^,:PlC.2i?P9.?,lr7.1,2F6.2, 1C7.?,4F6 .2» ?F7.i , 
llc* • 3, 1 F6. ?» lc  7. 3 ) 

1010  F07HST(1MH2A4,33X443AG"I4I 

1090  F0941T  C10F7.3,?X,1M‘) 

1091  F0RW-*T  (10F7.3,2X,1HA) 

IF  CI900Y.  tO.  D GO  TO  • 21 
9P*Y*  (2. *AX 2*92(7)  -?.*AXX*D? (4)  ) 

AP«?.*AXY*\/*  V-02(  3)  *91  (12) 

XI*A»CTAN(9P, AP) 

BETA--ALPHAN*SIN( XI) 

alpw=  alphan*cos(  xn 

50  T 9 A 4 1 

021  ALPH=ATAN(91  (14) /Old  2)»  •oosn 
9tTAr  AT  AM  (01  (1  3) /Ol<  12))*C0»I9 
941  «*p«3p*C0N9 

IF (“99 (LN, 50) .NE*  0)  GO  TO  2 
PGsPSM 

PRINT  101C»TITLt,P5 
PRINT  1000 
PRI N"  1001 
PRI»r.  1002 
2 p*  9! (5) 

XFCP.rO.O.)  GO  n 4762 

SS9C= (ni(5)*01C5)*A(15)*A(15))  / 1 0*5* A ( 19) • CHA* A ( 16 > ♦ 4. ) 

TAJ  * 1 • /SORT ( 1 • - 1./5S9C) 

029s  A(|9)*A(19)*Ad6) 

ClAKP  * 0*SI/C?.*  /)  € 1 p ) - V)  MDM*  M .»TAU)«  92.1/C  ?.  • A ( 16)  ) • (CF  (2 7 ) 

1M1.-7AU))  - 02M/ A(15)  •^‘)»a«TA!ll 
CL A KM  s (Q*S)/(?.* A ( 1 5) »V) M77AM t.-TAU) ♦ 3 2K/( 7 . • A ( 16 ) )• (Ct(?7) 


00004940 
00004950 
0OQC4S6Q 
00004970 
00004930 
00004690 
00 0 C 490  0 
00004910 
00004970 
00004930 
00004940 
00004950 
00804960 
QQG34Q70 
00004931 
0QCQ4V4Q 
00006003 
0CC  060 10 
00005320 
00006070 
00005040 
Q00Q5Q50 

00005060 
P0  0 06  o 7 0 
00005030 
0000r040 
00006100 


BEST  AVAILABLE  COPY 
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JUSWIMC  EJ  76/76  o«>r*t 


n»  <>.5*1.11. 


0V17/77  1T. 09.0C 


60 


*5 


7« 


75 


as 


65 


0M1.«TAU))  ♦ 071/A<15)*CMPA*T4'J> 

*IA1P«  IP*A<16()/<?.*A(1SIIM1.-  1./T4UI 

aa*  rtaMH'CiA-p-rp-wLHi'ji  •j’-vaaip) 

ULIHNs  <P«Afl6) 1/ <?.*» <151 (• ( l.»  l./TSU> 

69*-ClANP»  (P-WLH'II  -rL4’*'l,(>-WL  410) 

»53«  n»S*»C19»,:P<?C)/<«<t5»*5QET(**»»2  * <n,,7m»C0NO 
SO  TO  (.763 
6767  Cll«*l). 

tiling. 

HL  AHH  = C • 

»<3«3. 

taj.c. 

56900.0 
6763  CONTIN'IE 

“PINT  130  3,  <01  ( I)  ,I»1,4)  ,7,91  (5)  ,«L»H4M,»1C‘<N,',P,»LPH,9ET»1 
lCl»“M,!r6«1P,HL»11,iL»H»,S5>):,T4U,«K3 
IFIIOUN.FO.O)  GO  TO  1113 
LC0I3NT  * LCOON7  * 1 
IFaOOONT.EO.lO)  GO  TO  1117 
9»l(Lr0UNT)  =f.LPH 
06E 1L  G OUNT) = 9E  TA 
SO  TO  1113 

111?  0»ll<.r0UNT)*ALOH 
09E ILGOUNT) = 9ET  A 
LCOJNTsO 
1H3  insimi 
»ETU»N 
END 


60005110 
000  05 1 TO 
00905130 
000  0'  1<*0 
60005  1S0 
00006160 
00005170 
00005130 
60005190 
60006700 
00  C 0571 0 
00005770 
60005730 

00006750 
00905760 
00005770 
00005790 
00005790 
00005300 
00005310 
00005370 
000C5330 
00006360 
0P 1 05950 
0 05.350 
00006390 
000C5400 
00006610 


best  available  copv 


AFIT/CA/AA-77D-8 


SWCiriNE  IHIT6 


74/74  OPT *1 


pin  4.5»4i4 


SUBROUTINE  INI  TS 


C sitt?f»  ICO  C.  V • INGRAM,  AFRO-S®ACE  FNGR. 

C UNI  VERS  I T Y OF  NOTR?  *)«ME,  NOIRE  OAME,  INO. 

c 

C iNimiTMTION  PROGRAM  FOR  N-OFGREE  3F  FREEDOM  NOTION  C3»N*6) 

C 

REAL  -aSS, IX, IY,L 
tOOt  FORM*  T ( 1244) 

1002  FORMAT  C////1X#?^M£;H0  Oc  INITIALIZATION  INPUT*//! 

1003  FO^HST  (4(6<1X,?HA(  , 12, 4H)  s ,Fll.4)/)> 

1004  FOR*aT  <iXf7HNOR:  = , Il,3X,7wt£>RM  = , I2,3X,7HNALL  * .11,  3X, 
S/HIPIN  « ,11, 3X,  7HS3DOY  = , 1 1, 3X , 9HI SCALE  * ,11, 3X,6HI TRST  = , 

*/> 

1005  FORMAT  <//5<6< IX, 6MSCALF< ,12,4-D  r ,E6.2)/>> 

10  06  FORMAT  (/,5X,4HTlME,5X,6MTMRJSTfTX,4MH  ASS,  5X,UMCG  P0SITI0N.7X, 
1-02,2MIX,7X,2HIV,  /) 

10  0/  F0RM8TC4X,e5.2  ,6<  , -5 . 0 , 6X  , . 3 , 6 < ,F8  . ?,  OX,  F6.J  t 4XF6  . 3 , *,X  ,F<>.  3 > 

1006  FORHflT(5t,  13H7HP-JST  CATA  *,JX,6HEPS  = ,F5 . 2 , 3X  , 9HTMETAC  = ,F5. 
t,5METfl  S ,F5.2,3<  ,*HL  = ,F5.2,//> 

1009  F0?HaT(//,5X, 14H30DY  LENGTN  * ,F*.2,/) 

1010  rORM4  T ! 71 1 0 ) 

1011  cOR-ftT( 1 10,3  21 0.4,110,210.4) 

1012  FORMAT (9210.4  / 3210.4  / 6E10.4)  % 

1013  cORMA  T ( 6E1 0. 4) 

1014  FOR-A* (2210.4) 

1015  FORMAT! 3(3210.4  /) , 6E19.4) 

LOGICAL  LSCALE,I4LL,LTRST 

COMMON  /UIN9LR/  NOW,  YVH12S  ),  *(123,3) 

COMMON  /AIRBIK/  TEN, PRES, RH3,<VIS,VA 

COMMON  /0EN9LK/  NOO, Y O <123 ) , "FN < 1 23 ) , TEMP < 128) 


09/12///  17.09.0 


00000610 
0000*1620 
00000630 
00C0C640 
00091650 
00000560 
00000679 
00000690 
00000690 
00000/00 
00000/10 
00000/20 
It,/ 00 000730 
00000740 
00 C 09  75 0 
5HCG0C000763 
000007*0 
03C39763 
2,3X00002790 
00000800 
00600810 


00000620 

00C03829 

60000640 

00000650 


COMMON/ C OF  ° 1 / NpL  V 

COMMON  /POL/  CE  ( 3 0 ) , ALPHAR,CMAtCMt»ft,C*a 
COMMON  /IN09LK/  I OEN , I WINO, I T'R,  I PRN , I ALL , NORC 
1RST 

COMMON/I  NIT /TI TLE (12) , SCALE (33!  , A (22) 

COMMON/ THRUST/ TIM  <15 ) , T®ST ( 1 5) , MASS < 1 5> , CG < 15) 

S3) 

COMMON/TSTOaT/EPS, TMTC,ETA,L  ,9L 
REI0C5, 1001)  TITLE 

DO  10  1*1,15 
YZHCI) *0.0 
IRST ( I) *0  • 0 
MISS (I ) =0 .0 
CGCI) *0.0 
IV(I)«0.0 
IY<I)«0.0 
CGPCI 1*0.0 
10  CONTINUE 
I9EN=  0 
ININO*0 
OO  11  J*1 , 3 
OO  11  t* 1,128 
tfCIvJMO.O 
V«MI)*0.0 


60000660 

00000670 

,IPUN,N900Y,ITRST,NTC0033380 
000  0C890 
00099993 

,IX (15) ,IY (15),CGP(13QC32910 
00039920 
000009T0 
00030940 
00000950 
00003960 
00003970 
00009930 
00000990 
00001900 
00001010 
00001020 
60001030 
00001040 


00QQ1Q50 

00001060 

00001070 

00001030 


8®  W, 


mi  ®'1 


AFIT/GA/AA-77D-8 


subroutine  iNir»  rs/r 4 opt»i 


HH  4.SH.11,  05/12/77  17.05 


roiixo.o 

•0001099 

)'Ntn«o.o 

00001100 

to 

TEN»(I>«0.0 

00001110 

11  CONTINUE 

00001120 

e 

00001130 

C 4500  CONTROL  PARAHETER5 

00001140 

05 

REAS (5, 1010)  NORC,IPRN,NAlL,TpjN,N909Y,ISCALE»ITRST 
LSCALF=ISCALE.EQ.l 

00001160 

XALL*NALL • EQ.  1 

■0001170 

LYRST*ITRST • EO .0 

■eooitso 

C 

00001190 

70 

G 4E00  INITIAL  C0N0ITI9NS 

00001200 

IF(LTRST ) 00  TO  4 

00001210 

RfAOlS.lOll)  NTR5T,EPS,THTC,ETA,L,9l 

E®S=EPS/S7. 295779 

00001230 

TNTC*TmTC/57. 295779 

00001260 

75 

50  T9  7 

■0001250 

4 CONTINUE 

•0001250 

HTR$T«1 

10001270 

. 7 CONTINUE 

• 0001200 

00 

READ (5 ,1012)  (A(I), 1*1,22) 

IF(LSCALE)  GO  TO  1 

00001300 

DO  2 1*1,30 

■ 0001310 

7 5CALF(I>*1. 

•0001320 

50  T9  3 

10001330 

OS 

1 4E»0(5,1015)  <SCALE»I) ,1=1,30) 

5 CONTINUE 

10001350 

XF (LIDST)  CO  TO  5 

00001360 

45*9(5,1013)  TIH(l),“AS5(l),rO(l),IX(ll,IT(l),;GP(l) 

READ (5 ,1013)  TIN (13) ,NAS5< t9),C0(I9) ,IY(I0) ,XT(IO>,CG®(I3l 

40 

REAO (5,1014)  (TIN(I),TRST(I),I»1,NTRS)) 

5 CONTINUE 

•0001600 

- 

e 

00001610 

C ECHO  OUTPUT 

00001620 

PRINT  1002 

00001430 

45 

PRINT  1004,  NORC, IORN,  NA IL , I PUN, N900T , ISCALE.I TRST 

00001660 

PRINT  1003  , (I, A(  I), 1 = 1,22) 

00001650 

PRIN*  1005,  (I,5CAUE(1),X=1, 30) 

00001460 

IF(LTPRT)  GO  TO  5 

00001470 

EEPS=EPS*57. 295775 

■000)490 

100 

THTC=THTC‘57. 299779 

00001690 

PRINT  1005, EEPS.TT-tTC, ETA, L 

00001500 

PRINT  1009  ,9L 

00001510 

PRINT  1006 

00001520 

PRINT  1007,!TIM(l),TJSTU),NAS5(t>,CGP(I>,CG(I),IX(I),IY(I), 

I*t,NTC0001S30 

105 

JRST) 

0000)540 

6 CONTINUE 

00001550 

RETURN 

00001560 

SNO 

00001570 

AFIT/GA/AA-77D-8 


SlMftWTlMf  SIXOEG 


MN  4.5*414 


• V17/77  17.09. 


susr-iutinf  si*oc".Mriosi 

e university  9 r wr«  dahe,  notre  done,  ino. 

1001.  FORNATd/N  ALPHA  F9POR7XE14.9) 

isoo  fornatuim  wind  e»°d’> 

1006  eoRHATMin  HACH  E RR0R3XC14 . S ) 


0000169  0 
00001640 
00001609 
00001610 
00001670 

00001630 


1157  FO®H4T  (1X,3F70.5>  00001630 

0109  FOIUTI///  5x,» ALOHA  •>l«'10.6,»  OEGRCES  EXCEEDS  HAXIMUH*, 

1*  ll°HA  m A®  - AY •* 

1666  E0R4AT(//5X,. HSCH  ®,1F10.4,®  EXCEEDS  HAXIlUH  HACH  », 

9*NU“"FR  IN  ARRAY •> 

REAL  HASS, IX, IV,1TX, NTr,HTZ,L,HYJO,NZJO  00001640 

OEAt  "DOT  00001650 

INTESFR  PC  000016A0 

LOGICAL  ICZA.tCHt.LCHBA  00001830 

LOGICAL  LT»ST  _ 00991840 

COHHON  /AIROLX/  TEH, P®CS,®MD,XVIS,VA  00001670 

COHHON  /COEE/  AHCH(2C,  30) , AL°HA (7 0,3  0) ,C( ’0,70,301  ,49A (3  0)  ,494(30)  0000*600 
COHHON  /IN08L*/  I DEN, iHl 43, I T PR !, I PRN, I ALL , NOPC , 1 PUN.N900Y, I TRS T , NT  00 0 C 165 9 
IRST  00  0 C 170  0 

COHHON  /INIT / TITLE(17),S:iLF(!0) .A(Z?>  00001710 

COHHOH/THRuS  T/TIN(15),TRST(15).HAS$(15),C5(15),IX(15),IY(15) , CCP ( 1 0 3 0 C 1 770 
tl>  00001730 

COHHON/TSTOAT/EOS.THTC.ETI.L  ,71  00001740 

COHHON  /WIN4LK/  NOH, YH  (1701 , 4(129,3)  #0011750 

C04HON7COEP1/  4P. Y #0001760 

COH4DN  /POL/  CF < 3 0 ) , ALPHA®, rHA,I4PA,C?»  #0001770 

C0H4DH  /DEN?L</  NOO,  YD  (179)  , O-'N  (l  79)  , TEH»(  178)  00091790 

COHHON/ FD9X/D1 (1. I .0  7(141 , 971(14) ,DE?(14) .MINI  3) , DAL (10) , DOE  (10)  00001740 

C0HHPN/E09X7/V,D,  S.O»,Alphan,«hC-(n,»;,LN  00001900 

JOHHDN/E09X3/6XX, A3Y , A XZ  , AVX , AY Y , AY  7 , A7X , A7Y , A77  (0001910 

COHHON /E 09X4/  LCDUNT 

OATA  CONO.PAO/67. 295779, .017453749/  (0001870 

LT4ST«ITRST.EQ.O  00001850 

ItER=0  #0001860 

19*0  #0001870 

P5«0  (0001980 

KO’JNT*0  #0001099 

LCOUNT*#  00001990 

94*1  00001910 

07(11*1.  000019*0 

S*3.1415977*A(  19)  »»7/4.  (0001930 

10  70  1*1,14  (0001940 

70  01 (Z ) *A (I)  10001450 

RAJ?*RAO/7.  (0001960 

0*91 ( 8 ) *RA07  (0001470 

E*D1(9)*PAD7  00001990 

r*Ot < 1 01 *RA07  #0001490 

CP*COS(E»  #0007000 

SP*SIN(EI  (0007010 

CT*CDS(D)  00007070 

ST*SIN(0>  C (#007030 

CS*C06(F)  #0007040 

6S*SIN(F1  (0007050 

C VEPSDR  COHPONPNTS  (50)  (0007060 

31(81  *-CP*ST*SS*SP#CT*OS  100 07070 


BEST  AVAILABLE  COPV 


I 


92 


AFIT/GA/AA-77D-8 


SU4»*l*lNf  SIX0X3  76/76  0PT*1 


FTM  6.5*414  09/1 2/77  17.09.0 


oi(9»«-cp»ct*ss*5«>»5T*c5 
Dt(ie>*CP*sT*cs*SP*CT*$5 
M 9l«in*CP*CT*CS»SP*ST*SS 

I9P*A  (10) 

1 « A(20> 

25  lrmzoiMDim-*  (2in  .gt.o.igd  ro  597 

91*1 

•5  55  YOKJ) 

• IMITRST)  CO  To  56 

1*01(1) 

; IF«X.CT.TlH(MTRSm  50  TO  55 

I 1*1 

70  50  XFfX.LT.TlHClMM  SO  TO  57 

SO  TO  58 

‘ 57  20HTINMC 

! Y|«TR$T (I) 

• 7f  Y2»T7ST(I*l) 

▼■Y1MX-TIKI)  >•(  (Y2-Y1I/(TIM(IM)-TIN(I))  I 
‘ ATRST* Y 


00007080 
00007C90 
00007100 
00307110 
00007170 
00007130 
• 00*»21fc0 
00007150 
00007160 
00007170 
00007150 
00007190 
00007700 
00007710 
00002770 
00007730 
,00007740 
00007750 
00002760 
00007770 


Yl«H*SS(tl 

00  Y?*HASS  C9TRST) 

Y*Y1 ♦ (X-TIH(I) ) • ( (Y2-Y1)/(TIM(NT*ST)  -TIM(l) ) ) 
M16)*T 

9D0T«CAHP-4(18) )/A(20) 

8MP«a (18) 

•5  Y1*S6 (1 ) 

?2*CG (MTRST) 

?«Y1MX-TXH(1>1*  ( (Y2-Yl)/(TIH(NTRST)-TIM(l>>  ) 

ACS*Y 

Yi*lX(l) 

90  f 2*  t<  (MTP.ST  ) 

t«Vl* CX-TI9C1) )*( (Y7-Y1)/(TIH(NMST)  -TIK(t) ) ) 
I(151*Y 
»1*IY(1> 
f ?*IY (MTRST) 

95  Y*Y1MX-TIM(1)  )•( (Y7-Y1)7(TH(NTRST)  -TIM(i) ) ) 

0(16) *Y 
Y1*CCP(1) 

Y?*CCP( NTRST) 

Y*Y1MX-TIH(1) )•( (Y7-Yl)/(TIH(NTRST)  -TIH(t))) 

100  2GL*Y 

56  COHTTMMC 

ir(X.CT.TlM(NTR5T  ) ) *lOOT*0. 

5CL*1  ./S09T  (01  (8)  ♦•?♦  91  (9)  ••  ?♦'>!(  10)  ••201  (11)**7) 
Oll6)*5CL*r)l(8) 

105  01(9)*SCL*91(9) 

01 <10)*SCl*Ol (10) 

St(ll)=SCL*91(ll) 

CALL  ATMOS(Y) 

C RATRIX  Or  OIRICTION  20SINES 
110  AXX*01 (6)##7^0l(l l)*»2-.5 


000  0275  0 
00007790 
00007330 
00007310 
00007370 
00007330 
00007340 
00007350 
00002360 
00007370 
00007330 
00002390 
00002400 
00007410 
00C07420 
90007430 
00007440 
00007450 
00002460 
00007470 
00007430 
00007490 
00007500 
00002510 
00002570 
00007530 
00002540 
00007550 
00C07560 
00002570 
00007530 
00007590 


AXY«’31(8)*f)l(9)Oltl0)*ni(ll)  0000  760  0 
0X701  (8)*01C10)-01(1 1)»0!  (3)  00007610 
AY<*01 (ft)«Dl(9)-91(ll)*Ol(10)  00007620 


AFIT/GA/AA-77D-8 


SUWOl’tNF  StxOf  G 


7N/7*  0»*T»1 


FIN  *.$*<■!« 


•.717777  17.09.0 


115 

IYY«91(9>**2»0i(l 1) 

0000  26  TO 

* 

M?*01(9»*01tl  0)  *01!  1 t)*9l(4) 

•0002640 

f 

0*I*P1(S)  *9t(  I'M  *01  ( * l)»f)l  (9) 

00007660 

1 

o>*ni t io) -ot  ( ii>*3i(5) 

00002660 

477*"! C 10» •• 2»01 C 11) ••?-.$ 

00002670 

in 

irtiwim.EO.O)  GD  TO  2104 

00002690 

c m i no 

00002690 

1my.lt. ywu>)  ro  to  ioi 

, 

00007700 

: 

!•? 

00002710 

{ 

100  Ifttmn.GE.T)  CO  TO  210 

00002720 

1 . 1» 

1-1*1 

00002730 

) 

CO  TO  100 

00002740 

1 

til  PRINT  1005 

00002750 

1 

RETURN 

00002760 

l 

*10  *t$r*  (r-YM<i-i)>  /<r«<  D-Hd-tn 

000027T0 

t ISO 

OO  2’0  <*1,3 

'•0002790 

» 

220  MI <(«»**(!-!,<)  ♦*4Sf • -mi»i|<)  ) 

OOOC27RO 

1 

1 

C MIND  COMPONENTS 

00002003 

rfX*2.M  WlN(i>»  AX  ONI  NT?)  • A X Y ♦ Wf N ( 3) • A XZ) 

00C02510 

s 

WY*2.»(WlN(t>*AYX*WlN(  ,)*AYY*WTN( 3 ) * A Y Z) 

00002623 

i 1S5 

WZ«2.* (WIN < 1) • 4Z < ♦WI N( 3) *47Y ♦«!  <( 3)  * A Z 7) 

00002630 

j 

C MEtOCIT Y VECTOR 

00002640 

! 

01(1?>*01< 12)-WX 

00002950 

Dl(13)*Dl (t  3)-WY 

00002860 

01  (1  Jr  1 * 01  ( 14)  -WZ 

00002679 

1st 

C SPEED 

OCO  0 239  0 

2104  V* SORT  (01  f t2)**2*01(13>*»201(14>*»?> 

00002690 

r 

C ALPHA  PRIME  - YAW  - 4RCTAN(SQRM  V*V#W*W) /U>) 

00002900 

j 

AL*MA*«ARCTAN(S9RT(D1 ( 13) ••2*01 (14) •• 2) ,01 ( 12)1 

00002919 

iss 

alphan=alpmar*cono 

00002920 

C DYNAMIC  PRESSURE 

000C2933 

I 

}r#$»9HO*V*V 

90002949 

C CROSS-SPIN 

00002950 

0PP*S0RT<91  (6)  **2  01  (7)**2» 

' 00002960 

151 

C ROLL  ANCLES 

•0002973 

OP* APCT AN (01 (13).  91(14)) 

00002963 

OP** A PC  T AN ( -01 (7)  ,01 (6)) 

00002990 

C PO/2V  ANO  90/2 V 

00003090 

* 

IMV.GT.C.)  CO  TO  22 14 

0C00301C 

155 

2 CONTINUE 

00003020 

i 

- PO^O. 

000  0 TO  30 

00*0. 

•000 '340 

OMCMNrQ. 

• 00  0 3050 

50  TO  2215 

000C JC60 

160 

2214  P0M01  (5)*A(19) ) / (2.*  V) 

•000307Q 

09* (npD*A ( 19) )/(?.*V) 

•0093090 

IMCHNsV/VA 

00CS3C93 

22 15  0S0*O* S* A (19) 

•000  1100 

IS  Ns (9*  3)* 1 • / A (13) 

00003110 

165 

c interpolate  coefficient  tables 

•COOTI’O 

. 

00  lr0  <*1,33 

0000  3140 

IPMNOA  (K)  .EQ.O)  . ANO,  (NOM(<»  .E9.SU  GOTO  145 

10003169 

1*1 

• 000  1160 

J*  1 

I0C0317* 

in 

I F (NOA ( < ) • EO . C ) 50  To  1?0 

000  0 3190 

IFfALPNAN.LT.ALPMAd.O)  SO  T9  111 

00003190 

BEST 


94 


AFIT/GA/ AA-77D-8 


surrqtine  stxors  74/74  opt»i  ft*  4«5*6i4  §9/17/77  17.09. 


M*MOA  «> 

I*(AlPMAN.GT  . AL°^A  <N,  K) ) GO  TO  0888 
i«? 

17$  HO  IFfALPHA  (J.KJ.GE.  ALPHAS)  GO  TO  115 

J»JM 
SO  TO  110 

111  PRIN*  1004, alphas 
RFTURN 

181  89  88  PRIN*  8899,  ALPHAS 

RETURN 

11$  FRA*<  A\.Pv4A‘»-ALPMA  t J-l  ,X>»  /l  AL°-*Af  J,<>  - At »-« A ( J-l  , <>  ) 
170  IFfNOMfO .EQ.0>  SO  TO  140 

IF(A-CMN.LT.A‘1CH(  1#<)  ) GO  TO  til 
185  M-NO*  t K| 


10009700 

' 00009710 

00005770 
00001730 
00001740 
00001750 
00003760 


0000*770 

00003780 

00003790 

00001300 


IFfAMCMN.GT.AMCH<Nt*) ) GO  TO  1131 

1*7  60003310 

130  XFfAMCMfl.K) .GE.AHCHN)  GO  TO  135  00003370 

190  !*!♦ 1 00003330 


195 


700 


705 


710 


715 


770 


775 


GO  TO  130 

131  PRIST  10C6,AMCHN 
RETURN 

1131  PRIST  1665*  A*lCH*f 
RETURN 

135  FRMr  CAMCHS-AMCH(t-l , K) ) / C A ( T , < > -A MCM ( I -1 , <) ) 
:ECO*C<I-l,JtK>  *FRR*  (C(I.  j,r»-c«  X-l,  j,kh 
ircsOACO  .ET.0)  GO  TO  150 

vi*;u-i»  j-t.oppRip  <cu,  j-i,o-:ci-t,  j-i,<n 

LC?  A- K. E9.6 
ICNA*K.E0. 13 
LC^PA  *< • EO. 19 

IFimZA  ) C7A=  (CE  to  -TD/CALPKA  ( J*K)-ALPH4  < J-l , KM  *C0ND* 
SSCALf  OO 

IFCLOSA  ) CSAs(CF(<)  -T1)/(ALPMA(J*K) -ALP-<A  (J-l*<>  )*CONO* 
SSCALPtO 

tFCLCMPA)  C*1PA=  (3F(<)  -Tl)  /<  AL*HA  ( J#  <)  -ALPHA  (J-l,  <)  > *CONO* 
f 5CALF  CO 

SFOO  *T  1*FRA*  f CF  f K>  -Tl> 

SO  TO  150 

140  CEOO*C(I,  J- 1,0  ♦FRA*  { C (I , J,  O -C  < I , J - 1 ,K  >) 

SO  TO  150 
145  3E«K> *C(1,1,K) 

150  CFfO=*E(<)*SCALF  (<) 

C SOMPUTF  rorFFICIENTS  CX,  CT,  C7,  CL,  CM,  CM  (TA9LE  3) 

1511  SOSTTN'JE 

IFfNPir.NF.ll  GO  TO  1517 
OO  1^15  K<=1,10 
1513  3E(<0*CF<K0  *SCUFC<0 
1517  CO*!TI»»ut 

SP«  $IN  fOP» 

CP*C05C0P) 

SP4*«itN(4.*0P) 

CP4«00S  I4.*0PJ 
$®P4*SIN«4.*0PP) 

:pp4*cos(4.#oppj 

5x*cr m 


00003340 
000  0 3350 
00003369 


000C3370 
00001390 
60001390 
00003400 
C000.1410 
00  0 C 3470 
00001410 
00003449 
0000  *450 
0000  *469 
00003470 
0000*450 
00  0 0 349  0 
0C  0 0 *50  0 
8000  3510 
80003570 
00093510 
00  0 0 *540 
0000  1550 
00001560 
00003570 
09001580 
000  01590 
00003600 
00001610 
0000*673 
0009*630 
800  0 164 C 
00001650 
09003560 
00003670 
00003690 


AFIT/GA/AA-??D-8 


sust'xrm  sixoed  74/74  opt.i 


FTN  4. 5^414  09/1  2/77  IT.  09, 


no 


035 


040 


043 


050 


055 


050 


CYMCFf  2)  ♦CE(3)*SP«*CE(4)  •C®4^*0*CE(5) )•:»♦  (CC  < & > ♦ CF(7)*S®4  0E(3  )•  00  CC  *69  0 
fCP4)»S®^CEt?2)  COOOVCO 

5?*(Cf  (2)»CE(3)*5p**CE(4)*CP4*®VCE(5I)M-SPW  (CE  (61  ♦ CEC 7)  *S®4  ♦ CF  ( 00 C : *7i  3 
ld)*CP4)*C®»CE (23)  000Q?7?0 

CL*Cr  (9)  ♦ CEdO)  OE(  19)  -SOw^E  (It)  •C®4^®0»CE(12)  ♦ALPHA9*(CE(241  ♦ CE(CCC2  37  30 
105)®?P»CE  ( 26 ) •C®)  C0033743 


CM*  (CF(  13)  *CF  ( 14)  *SP44CF(  15)  •C®4»  *C®*  (CE  C 1 6)  ♦Cf  ( 17)*SP4^pc(1S)*CP4C9C3‘t75J 

i*pd*:fc  i9>  )*sp*cr<?7>4rEc?3)  •*;=>>«♦  of  (?9)  •cpp*5*cd*cos  n®® » ♦r.E (?e>  ?or:  T7  6 3 

CM*  (CE(13)  ^CE(  14)  *S®4  4CE(  15)  •C®4>  • <-S®>  MCEd5*)  ♦CEd?)  •S©4^GE  (15)  *CQC9*770 


|C»4»®n*CE  (19) ) <CE( 27) *CE (23)  * S©®4 *CE ( 2 9 ) *C®P4 ) *QO*  <- 

*:ccoi) 

l«'llMlNn.EO.O)  CD  TO  151 
C 100  *IN9  ®AC< 

oici») *01(1?) ♦ wx 

9l(lM*0l(lJ> ♦MY 
9lC14)r0l(14) ♦ MZ 
C KIWEM  AT  Xr.  AL  DELATIONS Hi PS 
C • VELOCITY  (9) 

151  92(?)=2»M,)1(12)*A*X01(13)*AVX01(14)*A7X) 
9?(3>«?.Mni(l?)*A<Y»0l(13)*A7<Ol(14)*A?Y) 

92(4>*2.M01(  12) * Ax? ♦01(1 3) • A Y* O 1 ( 1 4) *AZ2) 

C • MOMENTS 

0?  (5) *1 */A (15 ) • (1ST* CL ) 

0?(6)=1./A  (16)  *nST*CMOl  (5)  *Dl  (7)*(A  (16)-A  (15)  ) ) 

92(7) «1. /A (16) •(1SD*CM-01(5I •D1(5)*(A(16)-A(15)> ) 

C • VE*SO®  (10) 

02(d=-.5*  (-01(5)  *31  ( tl)-91(7)*9l  (9)  ♦OKS)  *01(10)) 
OOC9)=-*5*  (-01  (6)  *">l  ( 1 1)  ♦ Dt  ( * ) • 01  (* ) -91(5) *91(10)  ) 

92(10 >*-.5* (-9 1(7) *01(11) -Dl (6) *31(5) ♦OK  5) *91 (S>> 

90(11) *-.5* (01 (5) *91 (S>*3KS)*9i(9>*0l(  7) *01(10 ) ) 
ir(HOP.C.ME.O)  GO  TD  233 
6*0.*  U (17) -.307533686-5*91 (3)) 


SI N ( OP®) ) ♦GCfO  *759 
0030  *790 
00C03530 
CO  0 0 351  0 
0000*520 
00303530 
000 0 35 4 9 
00003850 
000  3 7860 
000  0 397  0 
00003380 
CO 0 0 3590 
00003990 
COO  0 3910 
0000  3920 
00003930 
00003940 
0OC0395O 
00903960 
00003970 
0C00  3950 
00003990 
000  0400  0 


59  TO  ?50 

230  G*?«*  A ( 1 7) 

265  C • FD»Cf 

250  0?(l‘’>=Q$M*CX-O*ATY-(9l(6>*91(l4)-0l(7)*9t(13)> 
92( t 3) * OS*1*CY-G*AYY-  (91(7)  • 9 1 (l  2)  -01  ( 5)  • D1  ( 14)  ) 
0?(14)*0SM*C7-GM7Y-  (01(5)  *01(1  3)  -91  (6)* 01  (12)) 
IF(ITOST)  GO  TO  251 

270  TX=ATPST*COS(EPS) 

?V*AT©ST*SIN(PP$) •SrM(TMTC) 

Y2*-ATPST*SIN (E©3 )*COS(TMTC) 

9TX*FTA*TX 
NTY*-(9L-CGL)*T7 
275  MT7*-(,*L-CGL)*TY 

NY J9*M00T*  01  (6)*‘»  **2 
M?JOsMOOT*9l (7)*L**2 
02(5) *02(5) ♦MTX 

92(6)  *02(6)  ♦ <MTY^MYJ9»0*S*C7*A9G)  / A(16) 

280  02(7)=02(7)  MMT7fM*j*)*Q«S*CYMCG>  /A(16) 

02d?l*92(12)*TX/ A(lft) 

9?(13)*02( 13) ♦ TY/A(1M 
D?(14) *02( 14) ♦T7/A (1 8) 

251  CONTINUE 

265  IF(LM.EO.O)  CALL  E9 


00004010 
00004020 
00004030 
00004040 
00004050 
00004060 
00004070 
OOCQ4CSO 
00004090 
00004100 
00004110 
00004120 
00004130 
00004140 
00004150 
000041S0 
000041 70 
00004190 
00004190 
00C04200 
00C04210 
0090422C 
00004230 
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SUBROUTINE  SIXOFG  74/74  0“T*l 


riM  4.s»4i4 


09/12/77  17.09.0 


C 9UNGE-KUTTA  CONSTANTS 
90  298  1=1,14 

GO  TO  (294, 295, 236, 297, 70), 91 
294  JEim*01C!> 

3E2(I)*02(I>*M 
E1*DE2(I)» .5 
SO  TO  298 
299  «1*0?(I)*.5=H 

5E2(I)*0E2CI)»4.»71 
SO  TO  298 

296  E1*0?(I)*H 

JE2m  = 0E2(I)  *E1»E1 
50  TO  298 

297  El«(0?2(I)*02(I)»HI/5. 

298  01(I)*9E1(I)»E1 
9l*Nl»l 

GO  TO  55 

70  SO  TO  (540, 550), 44 
940  <OUNT<KCUNT<1 

17(01 (*) -1 (22) ) 550,597,580 
950  I7(A4S(91(3)-A(22)).LT..01)  GO  TO  597 
9*-(91 (3>-A<22)>/02<3> 

94*2 

91*1 

C NEXT  TIN*  LINE 
SO  TO  55 

180  17  (M00{  KOUNT  » TPRN  ) *NE  • 0)  GO  TO  25 
CALL  EO 
SO  TO  25 
597  CALL  F9 
RETURN 
EN9 


00004740 

00004250 

00004760 

00004270 

00004280 

00004298 

00004300 

00004310 

00004320 

00004130 

00034340 

00004350 

00004360 

00004370 

00004330 

00004393 

00004400 

00004410 

00004420 

00004430 

00004440 

0000445(1 

00004460 

00004470 

00004430 

00004490 

00004500 

00004510 

00004520 

00004530 

00004540 

00004550 
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iunctioh  hcton  t 4/74  opt«i 


rTN  4.9*414 


•9/11/77  i 


t 


f 


11 


u 


II 


FUNCTION  «*CTAN(»i,PI) 
1,1,1 
1 171*1)  3>S.4 
> »*0TJN=4. 7113990 
9CTU9N 

4 «9CT»4«1.  5707963 
9ETU»N 

I 171*1)  5,6,5 
6 IF1PI.GT.0 .)  GO  TO  6 

»9CT«N«3. 161592/ 

9CTU»N 
• »9CT»N«,0 
9ET09N 

5 «4C*lT*NtPt/*I) 
tF<P?)  JO, 9. 9 

II  09CT9 4*190*3,1415917 
9ETU9M 

9 1711*0)  It, II, II 
11  190714*190*6. 1031553 
9ETU»N 

II  19CT1N*A9C 
9STU9N 

mo 


ncooioo 

UC00010 
1900  010 
A9CC0030 
I9C00040 
A9C09C50 
A9C300S0 
A9C00070 
MC00090 

A9CJ0090 
»»c oo id o 
HC00110 
UC  001*0 
A9C90130 
A9CQQ140 
190013159 
A9C00160 
A9C0017  0 
A9C00190 
MC00190 
A9C00I00 
A9C00I10 
A9COOIIO 


7.19. 
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